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NOTATION 


body cross-sectionai area 
body base area (at jc = C) 
planform area 

reference area (taken as Ajj for the comparisons of computed with experimental 
results) 

surface wetted area 

exposed wing planform area (2 panels) 
speed of sound 

semimajor and semiminor axes of elliptic cross section 


axial-force coefficient. 


Roo^r 


crossflow drag coefficient of circular cylinder section, 


Qn^^^cy^cy 


drag coefficient, 


lift coefficient, 


lift 

QoeAf 


pitching-moment coefficient about station at from nose, 

Pn 

normal-force coefficient, 7— 

RooAf 

local normal-force coefficient per unit length 

pressure coefficient, 

Qco 

f 

side-force coefficient, y 

body cross-section diameter 
cross force per unit length along body length 
potential cross force per unit length along body length 
viscous cross force per unit length along body length 


pitching moment 

QooAfX 


V 





axial, normal, and side force 
corner rounding for body cross section, — 
longitudinal and transverse apparent mass coefficients 
body length 

body aftersection length 
body nose length 

Mach number component normal to body axis, sin a 

free-stream Mach number 

pressure 

free-stream static pressure 

dynamic pressure component normal to body axis, sin^ a 
free^tream dynamic pressure, 
body cross-section radius or comer radius 
arc radius of ogive 
displacement of crossflow 
semispan 
time 

pV^ 

free-stream Reynolds number, — - — 

Reynolds number component normal to body axis. Re 
body volume 

velocity component normal to body axis, sin a 
free-stream velocity 
body width 
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X 




X reference length (generally taken as d for the comparisons of computed with 

experimental results) 

X axial distance from body nose 

axial distance from body nose to aerodynamic normal-force center (center of 
pressure) 

x^ axial distance from body nose to centroid of body planform area 

Xf^ axial distance from body nose to pitching-moment reference center 

Xgf axial distance from body nose to aerodynamic side-force center 

a angle of attack 

P angle of sideslip 

6 wing planform semiapex angle 

1 } crossflow drag proportionality factor (ratio of crossflow drag for a finite-length 

cylinder to that for an infinite-length cylinder) 

H viscosity coefficient of air 

p density of air 

0 angle of roll about body longitudinal axis 

Subscripts 


b body base 

cy cylinder 

Newt Newtonian theory 

nose body nose 

o equivalent circular body or cross section 

SB slender-body theory 

stag stagnation 
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SUMMARY 

An engineering-type method is presented for computing normal-force and pitching- 
moment coefficients for slender bodies of circular and noncircular cross section alone and 
with lifting surfaces. In this method, a semiempirical term representing viscous-separation 
crossflow is added to a term representing potential-theory crossflow. In computing and 

for bodies alone, slender-body theory is used for the term representing the potential 
crossflow. For bodies with thin wings and tails, the linearized potential method of Nielsen, 
Kaattari, and Pitts, modified for high angles of attack, is used. 

For many bodies of revolution, computed aerodynamic characteristics are shown to 
agree with measured results for investigated free-stream Mach numbers from 0.6 to 2.9. The 
angles of attack extend from 0“ to 180® for 3foo = 2.9 and from 0® to 60® for = 0.6 
to 2.0. 

For several bodies of elliptic cross section, measured results are also predicted reason- 
ably well over the investigated Mach number range from 0.6 to 2.0 and at angles of attack 
from 0® to 60®. As for the bodies of revolution, the predictions are best for supersonic Mach 
numbers. 

For body-wing and body-wing-tail configurations with wings of aspect ratios 3 and 4, 
measured normal-force coefficients and centers are predicted reasonably well at the upper 
test Mach number of 2.0. However, with a decrease in Mach number to = 0.6, the 
agreement for rapidly deteriorates, although the normal-force centers remain in close 
agreement. 

For Afoo = 0.6, 0.9, and 2.0 and angles of attack of 10®, 20®, 30®, 40®, and 50®, 
vapor-screen and oil-flow pictures are shown for many body, body-wing, and body -wing-tail 
configurations. When separation and vortex patterns are asymmetric, undesirable side forces 
are measured for the models even at zero sideslip angle. 

These side forces can be significantly affected by changes in Mach number, nose 
fineness ratio, nose bluntness, and body cross section. Generally, the side-force coefficients 
decrease or vanish with the following: increase in Mach number, decrease in nose fineness 
ratio, change from sharp to blunt nose, and flattening ot body cress section (particularly the 
body nose). Additions of afterbody strakes, wings, or wings plus tail produce much smaller 
or no appreciable effects. 
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CHAPTER 1 
INTRODUCTION 

Over the last several years, high angle-of-attack aerodynamics has increased in 'mpor- 
tance because of the demand for gieater maneuverability of space shuttle vehicles, missiles, 
and military aircraft (both manned and remotely piloted). Until recently there has been a 
general lack of analytical methods and aerodynamic data suitable for use in the preliminary 
design of most advanced configurations for flight to high angles of attack over a wide range 
of Mach and Reynolds numbers. There has been, however, considerable research leading to 
the development of methods for predicting the static aerodynamic characteristics of simple 
shapes, primarily slender bodies of revolution. 

Prior to the work of Allen in 1949 (ref. 1), most analytical procedures for computing 
the aerodynamic characteristics of bodies were based on potential theory, and their useful- 
ness was limited to very low angles of attack. Allen proposed a method for predicting the 
static longitudinal forces and moments for bodies of revolution inclined to angles of attack 
considerably higher than those for which theories based only on potential-flow concepts are 
known to apply. In this method, a crossflow lift attributed to viscous-flow separation is 
added to the lift predicted by potential theory. This method has beeii used quite success- 
fully to compute the aerodynamic coefficients of inclined bodies (e.g., refs. 1 -6), although 
most data available for study until 1961 were for bodies at angles of attack below about 
20°, and the formulas were initially written to apply only over about this angle-of-attack 
range. 

In 1961, Allen’s concept was adapted by Jorgensen and Treon (ref. 7) for computing 
the normal-force, axial-force, and pitching-moment coefficients for a rocket booster 
throughout the angle-of-attack range from 0° to 180°. Reasonable agreement of theory with 
experiment was obtained for a test model of the rocket booster over the Mach number range 
from 0.6 to 4. The Allen concept was further applied by Saffell, Howard, and Brooks 
(ref. 8) in 1971 in a programmed method for predicting the static longitudinal aerodynamic 
characieristics of low aspect-ratio missiles operating at angles up to 1 80°. 

In 1958, a method for computing the aerodynamic characteristics for bodies of non- 
circular cross section was proposed by Jo';.ensen (ref. 6). In this method, normal-force and 
pitching-moment coefficients (C)y^ and computed by Allen’s formula*' for the 
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equivalent body of revolution which has the same axial distribution of cross>sectionaI area as 
the noncircular body. Then Ihe values of Cj^ and for the noncircular body are computed 
from and ratios determined from apparent mass coefficients (i.e., from 

slender-body theory). Good agreement of theory with experiment (ref. 6) was obtained by 
this procedure for bodies of elliptic cross section at the conditions investigated (alb - I 
to 2, 0 = 0" and 90°, M^=2to 4, and a = 0“ to 20“). 

Early in 1973, the Allen concept was again applied by Jorgensen (ref. 9) to further 
develop an engineering-type procedure for computing normal-**; rce, axial-force, and 
pitching-moment coefficients for slender bodies of circular and noncircular cross sections at 
a = 0® to 180®. The and formulas were written, however, for a body whose cross- 
sectional shape remains constant over the body length, but the cross-sectional area, of 
course, is allowed to vary. 

In 1973, Jorgensen (refs. 10 and 1 1) rewrote the and expressions to apply for 
the general case of a slender body alone or with lifting surfaces where the cross-sectional 
shape, as well as the cross-sectional area, is allowed to vary along the body length. For the 
si>ecial case of winged-elliptic cones, simplified expressions for and were also pre- 
sented. Good agreement between predicted and experimental results was shown. However, 
experimental results available for comparison with the method were limited to simple ellip- 
tic bodies and winged-elliptic cones at angles of attack less than about 20® and Mach 
numbers only from 2 to 4. Thus, it was concluded that additional comparisons for these and 
more diverse configurations were needed at higher angles of attack and lower Mach numbers 
to determine validity limits for the method. 

In 1972-73 there was a great need to enlarge the relatively small data base for bodies 
alone and in combination with wings and tails at subsonic, transoniv, and supersonic Mach 
numbers. This need still exists today, but it has been alleviated significantly by some recent 
introductory investigations into this higli a field (e.g., refs. 12-21). Most studies have been 
initiated primarily for application to missile aerodynamics. 1 wo recent studies, however, 
have been more aircraft oriented in that they have been directed toward the determination 
of experimental aerodynamic characteristics for slender bodies with thin wings (ref. 20) and 
with wings and a tail (ref. 21). There is still need for study and analysis of much of this 
recently obtained data both for slender bodies alone and for bodies with wings at very high 
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angles of attack. As previously mentioned, the Jorgensen method should be further com- 
pared with experiment. Also, there is an urgent need for contiiiucd study of unwanted side 
forces and yawing moments which have been shown to develop for models at high o valu«.*s 
with ^ = 0°. These unwanted side forces and yawing moments, which probably lead to 
undesirable stall/spin characteristics, have been m''-'sured on noses alone, bodies alone, 
bodies with wings, and bodic '. with wings and a tail (e.g., refs. 12-21). 

In view of the foregoing, the present study was initiated to accomplish the following 
objectives: 

1. Review and extend the derivation of Jorgensen’s engin.ering-type method 
(refs. 9-11) for computing the normal-force and pitching-moment characteristics of slender 
circular and noncircular bodies alone and with thin lifting su. faces. 

2. Assess the method for predicting the aerodynamic characteristics of missile-like 
bodies of revolution. 

3. Assess the method for predicting the aerodynamic characteristics of bodies of ellip- 
tic cross section. 

4. Assess the method for predicting the aerodynamic characteristics of bodies with 
thin wings alone and in combination with a tail. 

5. Present and discuss visual observations of the vortex flow fields over models of 
bodies alone and in combination with wings. 

6. Present and discuss data pertaining to the origin and alleviation of undesirable side 
for..:s and yawing moments associated with high a flight at ^ = 0°. 

In connection with objectivei (5) and (6), the present study is limited in scope to an 
experimental study of the vortex flow fields and side forces associated with bodies and 
body-wing combinations. It is firmly believed that a larger experimental foundation must be 
acquired at this time to aid in the further development of flow-field modeling techniques 
that can be applied with confidence in preliminary design studies. ThL< is not to say that 
important semiempirical modeling techniques hav not already be>’-. ‘Uudied. Most of these 
techniques have been based on the impulsive flow analogy and ar p ' • .lender bodies of 

revolution. In this analogy, it is assumed that the crossflow plane itwept uniformly down 
the length of an inclined body at the r'te cos a. Then theie is an analogy made between 
the leeward flow field and the developing wake behind an impulsively started cylinder in 
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crossflow. This analogy was initially suggested by Allen and Perkins (ref. 2) and applied by 
Kelly (ref. 22). 

The impulsive flow analogy probably has been used most elegantly in recent studies by 
Wardlaw (refs. 23 and 24). In one of his latest studies (ref. 24), he simulated the viscous 
crossflow plane by superimposing a large number of point vortices on the potential solution 
for flow about a cylinder. For some pointed-nose bodies, he has shown qualitative agree- 
ment with experimental observations, but even the latest Wardlaw method still must be 
considered to be in an early exploratory phase. Others have used the impulsive-flow analogy 
and attempted to model the crossflow field (e.g., refs. 25-32). 

Attempts also have been made to solve the crossflow field problem with the Navier- 
Stokes equations (refs. 33 and 34), but these studies have been applied only to bodies at 
supersonic speeds and moderate to low angles of attack. Computer times appear to become 
prohibitive for most pr;<ctical cases. 

Because of many deficiencies in the understanding of the physical flow fields around 
bodies and wings at very high angles of attack, further experimental investigation is required 
before much reliability can be realized from further analytic modeling. 
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CHAPTER 2 

DERIVATION OF BASIC METHOD FOR COMPUTING Cj^ AND 

CHARACTERISTICS 

Here we first review the derivation of some basic equations for computing the 
normal-force and pitching-moment coefficients, and C^, for slender bodies of revolu- 
tion to very high angles of attack (section 2.1). Then we extend the derivation to obtain 
more general Cjy and equations that can be used for slender bodies of circular and 
noncircular cross section alone and with thin lifting surfaces (section 2.2). Necessary empir- 
ical input values of crossflow drag coefficient for circular cylinders are presented in sec- 
tion 2.3, and theoretical formulas to obtain input local normal-force coefficients for non- 
circular configurations are presented in section 2.4. Finally, in section 2.5, we briefly discuss 
the relative influence of derived potential and viscous crossflow terms used in the basic 
equations to compute both Cyy and as a function of angle of attack. 

2.1 Bodies of Revolution 

In 1949, H. J. Allen (ref. 1) proposed a heuristic concept for predicting the static longi- 
tudinal forces and moments for bodies of revolution inclined to angles of attack consider- 
ably higher than those for which theories based only on potential-flow concepts are known 
to apply. In this concept, a crossflow lift attributed to viscous crossflow separation is added 
to the crossflow lift predicted by potential theory. For the potential lift, Allen used the 
slender-body equation derived in 1923 by Max Munk for airship hulls (ref. 35). 

From the momentum consideration, Munk (ref. 35) showed that, for slender (high- 
fineness-ratio) bodies, the potential cross force per unit hngth fp at any station along the 
body is given b'; 

fp=(K2- K,)q^smlx>L^ ( 2 . 1 ) 

where A is the cross-sectional area of the body at any axial distance x from the nose apex;ot 
is the angle of attack; and and K\ are, respectively, the transverse and longitudinal 
apparent mass coefficients for the body. From calculations of Ki and Kx made initially by 
H. Lamb (ref. 36) for ellipsoids of various fineness ratios, Munk (ref. 35) has shown that 
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K 2 ~ Ki is approximately unity for high fineness ratios. It has been customary to assume a 
value of unity for the factor K^- K^, and it is assumed in this formulation. 

G. N. Ward (ref. 37) has shown that the section cross force /p at small angles of attack 
acts as an angle midway between the normal to the axis of revolution of the body and the 
normal to the free-stream velocity (i.e., at an angle a/2). With this consideration, equa- 
tion (2.1 ) is multiplied by cos (a/2) to determine the normal-force distribution. At high a, the 
potential lift becomes small in comparison with the viscous crossflow lift, and it is of little 
practical consequence whether this low a multiplier [cos (a/2)] is retained or replaced with 
unity. 

Now consider the derivation of the crossflow lift attributed to viscous crossflow separa- 
tion. Allen (ref. 1), for a body of revolution of high fineness ratio, first treated each circular 
cross section as an element of an infinitely long circular cylinder of the same cross-sectional 
area. With this assumption, the local cross force per unit length due to viscosity is given by 

fy - ^ 2 . 2 ) 

where r is the body radius at any station x from the nose apex; is the velocity normal to 
the longitudinal axis; p is the mass density; and Cj is the steady-state crossflow (or normal) 
drag coefficient based on diameter, and unit length at stations. Cj is a fv ;tion of 
both the Mach number and Reynolds number components normal to the cylinder longitu- 
dinal axis. Hence, for a body at angle of attack, Cj is a function of 

“n 

= -J- = ^oo sin a ( 2 . 3 ) 

and 

Re^ = Re sin a (2.4) 

where 


Vn = sin a 


(2.5) 


Here ®nd are free-stream values, respectively, of Mach number, Reynolds 

number, and velocity; the speed of sound is denoted by a. 
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For this study, Cj is taken as the steady-state crossflow drag coefficient for a circular 
cylinder. Other researchers (e.g., refs. 22-32) have assumed that the development of the 
crossflow with distance along an inclined body of uniform diameter is analogous to the 
development with time of the flow on a cylinder impulsively set in motion from rest. In this 
analogy, it is assumed that the crossflow plane is swept uniformly down the length of an 
inclined body at the rate 1^ cos a. When this analogy is assumed, the value of Cj starts 
from zero (at zero time), then increases to a maximum value about 25 percent higher than 
the steady-state value, then decreases to the steady-state value for laminar flow (see, e.g., 
sketch (a) from ref. 27). For turbulent flow, experiments conducted by Sarpkaya (ref. 27) 
indicate that the steady-state condition is reached almost at the start of the motion. Thus, 
for turbulent flow both the steady-state assumption and the impulsive-flow assumption give 
the same result. However, for laminar flow there is a difference, but this difference is not 
studied in this report. 



Because of spillage flow around the ends of a finite length cylinder, the value of Cj is 
less than that for an infinitely long (truly two-dimensional) cylinder in the same free 
stream. Equation (2.2) should be multiplied by a proportionality factor q, which is the ratio 
of Cj for a finite length cylinder to that for an infinite length cylinder. This factor ij, 
which approaches unity as the cylinder length to diameter approaches infinity, is given from 
experimental results (discussed later). In the practical use of r\, it is assumed that the ijd of 
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the cylinder is the same as that for the body of revolution being considered. It is also 
assumed that the net-force effect of the front and rear end flows is approximately the same 
for both configurations. Thus, equation (2.2) is modified to give 


fv = sin^ a 


where 


1 


2 °° 2 sm^ a 




sin* a 


( 2 . 6 ) 


(2.7) 


Following Allen (ref. 1), the potential solution of Munk and the viscous solution are 
combined to determine the cross-force distribution along the body of revolution. With the 
potential cross force per unit length acting at a/2 from the normal to the free-si'^am 
direction and the viscous cross force per unit length acting normal to the longitudinal ax. . 
the total cross foice per unit length normal to the longitudinal axis is given by 


/ = /pCOS^ + /, 


or 


= sin 2a cos ^ + 2 tjC^ (sin* a)r 


a di4 


( 2 . 8 ) 


(2.9) 


From equation (2.9), the equations for the normal-force and pitching-moment coeffi- 
cients for bodies of revolution can be easily derived. With the normal force given by 
f dx and the normal-force coefficient defined by = (F^lqf^y), we obtain 


Ca = 


_ sin 2a cos (a/2) 




Likewise, with the pitching-moment coefficient defined by 

ri 

^ _ pitching moment _ •b ~ 

^ - -- ■ “ 


'm 






( 2 . 10 ) 


we obtain 


C. 


_ sin 2a cos (a/2) 
"» AyX 


r 


dA 


2t}C^ sin* a p 


r(x^ - x)dx 


( 2 . 11 ) 




where A^. is the reference area; X is the reference length; C is the body length; and is the 
axial distance from the nose to the pitching-moment reference center. 

General ii tegrated expressions for bodies of revolution at a = 0° to 180° can be 
obtained from equations (2. 1 0) and (2.11). For the sign convention in sketch (b), Jorgensen 
(ref. 9) has written the following equations for the normal-force and pitching-moment coef- 
flcients: 


sin 2a' cos ^ -I- tjCj si 

Af 2 “n Af 


sm* a' 


0° <a< 180° 


( 2 . 12 ) 


■^m 




0°<a<90° (2.13) 


and 


C = 
'-'m 




90° <a< 180° (2.14) 


where Aj^ is the body base area; Ap is the planform area; V is the body volume; is the 
distance from the nose apex to the centroid of planform area; and a' = a for 0°< o < 90° 
and a’ = 1 80°- a for 90° < o < 1 80°. 
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The axial distance from the nose apex to the aerodynamic force center is then given by 

X (2.15) 


'ac 


_ /^m \ 


When lift and drag coefficients are desired, they may be obtained at all a. values from 
the conversion expressions; 

Cl = Cjy cos a - sin a (2.16) 


Qj = Cjy sin a + cos a (2.17) 

Generally, it can be shown that for values of a well removed from a = 0° and 1 80® precise 
prediction of is not necessary to obtain reasonably accurate values of C^r and Cj^. 

Jorgensen (ref. 9) has suggested that, for rough engineering estimates, the axial-force 
coefficients for slender bodies can be approximated by 

Ca ; 0® < a < 90® (2.18) 

and 

Ca * cos* a ; 90® < a < 180® (2.19) 

Here cos* a is merely the ratio of the dynamic pressure in the axial direction to the dynamic 

pressure in the free-stream direction. A more precise method is also outlined in reference 9, 

along with procedures for determining Ca o and Ca . for often used conical and 

^a=o ^a = i80 

tangent ogive noses and flat bases. 


2.2 Bodies of Circular and Noncircular Cross Section 
Alone and With Lifting Surfaces 

For the general case of a body alone or with lifting surfaces where the cross-sectional 
shape can vary along the body length, procedures similar to those used for a body of 
revolution are assumed. Both the potential and viscous separation crossflow terms in equa- 
tions (2.10) and (2.1 1) for C/y and are generalized further. 

The potential crossflow term is generalized by multiplying the value inside the integral 

by the ratio iCJC„ ) - the ratio of the local nonnal-force coefficient per unit length C„ 

" "o SB r B n 

for the desired cross-sectional shape to the similar coefficient C„ for the equivalent circular 
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shape having the same cross-sectional area. The necessary ratios can be determined from 
apparent mass coefficients (slender-body theory) for many cross-sectional shapes. In 1958, 
this procedure was first shown (ref. 6) to have merit in predicting experimental Cj^ and 
results for slender bodies of various cross sections at supersonic speeds and angles of attack 
up to about 20°. However, because of very limited experimental data, further evaluation 
was not possible at that time. 

The viscous crossflow terms in equations (2.10) and (2.1 1) are further generalized by 

multiplying the values within the integrals by {CJC„ ) . Here the local ratio of C„ to 

" Newt " 

C„ at each x station is assumed to be given by Newtonian impact theory. Cj remains as 
"o “n 

the crossflow drag coefficient for the equivalent circular cylinder section. 

With equations (2.10) and (2.11) generalized as discussed, and for positive cbl/dA: 
values. 


C^- 


_ sin 2a cos (a/2) 


2t>C^^ sin* a 


and 


_ sin 2a cos (a/2) 
A^X 


sin* a 


/■&L " 

ra 
r (4L 


dx 


(2.20) 


(x^ - x)(bc 


r(x^ - x)dx 


( 2 . 21 ) 


In equations (2.20) and (2.21;. the first terms (from slender-body theory) are not applicable 
as written for winged-body sections where the dA/dx values are zero or negative, and 
procedures similar to those suggested in reference 38 probably should be used. Further 
adaptation of this method for use with body-wing and body-wing-tail configurations is 
considered in chapter S. 

In the second term of equations (2.20) and (2.2 1 ), there is some experimental justifica- 
tion for formulating the ratio iCjC„ ) from Newtonian theory and multiplying it by the 

‘ "o 

available experimental or theoretical crossflow drag coefficient for the equivalent circular 
cross section Cj . For subcritical crossflow Mach and Reynolds numbers, Jorgensen 


II 




(ref. 1 1) has shown that CJC„ values from Newtonian theory agree reasonably well with 

n 

those from two-dimensional tests (refs. 39-43) of elliptic cross sections and square cross 
sections with rounded comers. Jorgensen’s comparisons are shown in table 1 . Good agree- 
ment also can be expected at high supersonic and hypersonic crossflow Mach numbers 
where Newtonian theory by definition should be most applicable. The most doubtful 
regimes include the transonic crossflow Mach number regime and the supercritical Reynolds 
number regime. (These regimes are discussed further for circular cross sections in the follow- 
ing section.) 

Where reliable experimental crossflow drag data exist for a desired noncircular cross 

section, these data, of course, can be used. Then the values of Cj for the particular cross 

sections can be substituted in equations (2.20) and (2.21) in lieu of the product 

Ca {CJC„ ) , where Cj as now written is for a circular cross section only. Of course, 

S " Newt 

if the shape of the noncircular cross section varies along the body length, values of must 

“n 

be substituted within the integral in the second term of equations (2.20) and (2.21), and a 
great deal of empirical input data may be necessary for some configurations. 


TABLE 1.- Cj and CJC„ VALUES FOR TWO-DIMENSIONAL CYLINDERS OF 

" "o 

VARIOUS CROSS SECTIONS AT a = 90“ AS COMPUTED BY NEWTONIAN 
THEORIES AND MEASURED AT SUBCRITICAL MACH AND REYNOLDS 
NUMBERS 



NOTE : ALL C4„*l IN TABLE ABE BASED ON WIDTH OF 
CBOSS SECTION, NOT EQUIVALENT i. 
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2.3 Empirical Input Values 

2.3.1 Crossflow drag coefficient- To compute and for the equations derived in 
the previous section, values of crossflow drag coefficient Cj are needed for an “infinite 
length” or truly two-dimensional circular cylinder placed normal to an airstream. As pre- 
viously mentioned, Cj is a function of both the Mach number and Reynolds number 
components normal to the cylinder longitudinal axis, and hence for a configuration at angle 
of attack it is a function of sin a and Re^ = Re sin a. In these simple relations, 

introduced previously as equations (2.3) and (2.4), is commonly called the crossflow 
Mach number and Re^t the crossflow Reynolds number. For circular cylinders, necessary 
“state-of-the-knowledge” plots have been prepared for the variation of with M„ and 
Re„ (figs. 1-3). 

Figure 1 gives the variation of with Af„ over the range from 0 to 8. It was 

prepared from the data of references 44 through 49 and from data obtained recently by 

John M. Macha in the Ames 2- by 2-Foot Transonic Wind Tunnel. Also shown for reference 

are the theoretical variations predicted from Newtonian and modified Newtonian theories. 

Because of the close agreement of the Newtonian values with experiment at the higher Mach 

numbers, it is not surprising that computer programs utilizing Newtonian theories have been 

used successfully to predict space-shuttle-booster results in wind tunnels at hypersonic Mach 

numbers (see, e.g., ref. SO). Except for the transonic range, where data are very limited, the 

variation of Cj with M„ is well documented in figure 1 . 

“n " 

In the transonic range, the black symbols in figure I represent values of Cj obtained 
recently from pressure-distribution tests of circular cylinders of various diameters (1.9 to 
S.l cm) at crossflow Reynolds numbers from about 1.3X10^ to 4.9X10*. Many values of 
Cj were initially computed from the extensive pressure distributions measured by Macha 
on the cylinders in the Ames 2- by 2-Foot Transonic Wind Tunnel. Because there was a 
general increase in C^^ with decrease in cylinder diameter d (but not Re^) for values 
from about 0.9 to 1 .2, plots of Cj vs. d were constructed, and values of Cj were obtained 
by extrapolating the curves to d = 0. The black symbols in figure 1 represent these extrap- 
olated values of Cj^, values that should come closest to representing data for no interfer- 
ence between the models and the wind tunnel. These data agree well with the rocket 
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flight-test results obtained in 1953 by Welsh (ref. 49), but these or similar-sized models 

should be tested further in a larger transonic wind tunnel. 

As shown in flgure 1, there is a critical crossflow Reynolds number effect that can 

drastically lower the values of at below about 0.5. For less than about 0.5, if 

the crossflow Reynolds number Re^ exceeds the critical value of about 2X10* , the value of 

Cj decreases considerably. This variation is shown in greater detail in figures 2 and 3. 

“n 

Figure 2 gives the variation of with Re^ for less than about 0.4. It has been 
well documented over the last 60 years (e.g., refs. 39, 43, 45, 5 1 ) that Cj as 1 .2 for laminar 
boundary-layer flow and separation just before the critical Reynolds number of about 
Re^ = 2X10*. At about Re = 5X10* there is evidence (e.g., refs. 52-54) of laminar 
boundary-layer flow around the front of the cylinder to an angular position of about 80° or 
90°, where the flow separates, undergoes transition, and reattaches at an angular position of 
about 110° to form a laminar separation bubble. Then the turbulent flow separates at some 
position downstream (an angular location of about 130°). With a further increase in Rey- 
nolds number into the supercritical regime, the bubble decreases in size until the transition 
to turbulent flow moves upstream of the location of laminar separation, and the bubble 

disappears (ref. 54). From the low Cj value between about 0.15 and 0.30, Cj increases 

“m 

gradually, at least for an increase in Re^ up to about 5X10*. The supercritical Reynolds 

number regime has only been investigated recently in any detail (refs. 43 and 52-54), and 

there is still considerable uncertainty in the magnitude and trend of withRe^ andAf^. 

The shading in figure 2 indicates the approximate spread or uncertainty in Cj based on 

“n 

known data. 

Jones, Cincotta, and Walker (ref. 54) probably have made the most detailed study of 
circular cylinders in supercritical flow. With the use of freon gas to obtain high Re^, they 
have shown that there is an effect of on the variation of C^^ with Re^. Figure 3 (taken 
from ref. 54) summarizes their results for = 0.25 to 0.50. The reader is referred to 
reference 54 for their interpretation of these Cj results based on pressure-distribution and 
visual-flow studies. 

For noncircular bodies, as mentioned in section 2.2, experimental values of Cj for 

“n 

noncircular instead of circular cross sections can be used with slight modification to equa- 
tions (2.20) and (2.21). These data are generally not available for crossflow Mach numbers 
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above critical. However, some data are available for subcritical crossflow Mach numbers, and 

these data should be used if the crossflow Reynolds number exceeds the critical value 

(the Re„ where C,: drops rather drastically with a slight increase in Rc„). In table 2 some 
n Ufj 


TABLE 2.- REFERENCES FROM WHICH EXPERIMENTAL VALUES OF Q VS 

“/! 

Re„ CAN BE OBTAINED FOR VARIOUS CROSS SECTIONS AND FLOW 
DIRECTIONS 



original 

OP POOR QUAUW 
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references art l.bted from which experimental values of v.r jsRe^ can be obtained for 

variou' cross st'ctions and flow directions. It should be nt ,. most experimental values 

of are based oi, cruss-sectional width w and must tv - ij'-uplied by w/d, where d is the 
“m 

equivalent diameter ol fl’v' ctois section. 

2.3.2 CrossfloY) drag i»‘‘op >rti(>fi::lity jactf.- • ;i the equations to compute Cyy and 
Cf^, r) is the crossflow drag proportionality facior, that is, the ratio of the crossflow drag 
coefficient for a finite-length cylinder to that *ri ;»n infinite-length cylinder. Cylinder drag 
coefficients from which values of 17 can be determined have been measured (to the author’s 
knowledge) only at very low subsonic Mach numbers (refs. 56 and 57). 

Values of 17 for circular cylinders at very low crossflow Mach numbers (from ref. 56) are 
plotted as a function of length/diameter ratio in figure 4. Values of t} for flat plates are also 
plotted, but as a function of plate length/width ratio. The values for the flat plates are only 
slightly less than those for the circular cylinders. Thus, it is likely that figure 4 can be used 
to estimate values of rj for many cross sections varying from circular to flat. However, these 
values may be acceptable only for very low crossflow Mach numbers. 

An indication of the variation of 17 with crossflow Mach number can be obtained 
by computing values of 17 from high-o data (ref. 16) for slender bodies of revolution. 
From equation (2.12), 



Cjy - sin 2o cos (o/2) 
(ApMy)sin* e 


( 2 . 22 ) 


For two bodies of fineness ratio 10 and 12 (sketched in fig. 5), the variation of t 7 Cj^ 
with (for » 0.4 to 1 .6) has been computed from equation (2.22) with the use of 
data for values of a from about 45° to 60°. As shown in figure 5, the results for the two 
bodies agree closely. Now with the variation of C,i with M„ in figure 1 and the variation of 
^d„ ^>sure 5, the variation of 77 with has been computed and the results 

plotted in figure 6. (The circular symbols denote the computed values over the A/^ range 
from 0.4 to 1.6; the square and diamond symbol: represent values of 17 fo' vfry low 
cross-flow Mach numbers obtained from fig. 4.) 

For most supersonic and hypersonic values of Af^, q probably can be assumed to be 
unity, an assumption indicated as being essentially correct from past investigations (e.g., 
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refs. 3, 6, 7, and 9). TliC greatest uncertainty appears to be within the transonic range, 
and here further research is desirable, in this study, however, figures 5 and 6 are used in lieu 
of better information. 

2.4 Formulas and Values of (C„/C„ ) and (C„/C„ ) for Various Cross Sect! ns 

" "o SB <> Newt 

To use equations (2.20) and (2.21) for computing Cyv and C,^, ratios of the section 

normal-force coefFicients (C„) to those for the equivalent circular sections (C„ ) must be 

" "o 

determined. The ratios (CJC„ ) from slender-body theory are used in the first term of 

” "o SB 

each equation, and the ratios (C„/C„ ) from Newtonian theory are used ir ‘he second 

" "o Newt 

term. 

Formulas of (CJC„ ) and (CJC„ ) are now presented for some of the more 
o SB ^oNewt 

general cross sections encountered in missile and aircraft aerodynamics. Then, for several 

sample cross sections, values of (CJC„ ) and iCJC„ ) are plotted and compared. 

n .no SB ” "o Newt 

2.4.1 (CJC„ ) formulas- From slender-body theory (e.g., refs. 58-61, the ratio of 

n Ho gg 

for a winged-body cross section to that for the equivalent (same area) circular-body 

cross section can be determined for many cross-sectional shapes. In ref. II, (C„/C„ ) 

" "o SB 

expressions are determined for winged-circular and winged-elliptic cross sections (see 


sketches (c), (d), and (e)). 


[— S — I 
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Sketches (c). (d). and (e) 


For a winged-circular cross section with the wing planform perpendicular to the cross- 
flow velocit; (sketch (c». 




(2.23) 
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For a winged-elliptic cross section with the semimajor axis a and wing planform per- 
pendicular to the crossflow velocity (sketch (d)). 


where 





-I- a^) 


(2.24) 


ki — Oj - 


(a + bf 
4a I 


^ + h* ~ ) 


For a wingcd-elliptic cross section with the semiminor axis b and wing planform per- 
pendicular to the crossflow velocity (sketch (e)), 


where 



For an elliptic cross section without a wing (e.g., ref. 37), 



= X cos* ^ -I- - sin* <t> 
b a 


(2.25) 


(2.26) 


where ^ is the angle of roll about the body longitudinal axis, being 0“ with the semimajor 
axis a perpendicular to the crossflow velocity and 90° with the semiminor axis b perpendic- 
ular to the crossflow velocity (see sketches (d) and (e)). 


2.4.2 (CJCy. ) formulas- From Newtonian impact theory, {CJC„ ) 

” "ojVewr « ^o^■ewt 

expressions also have been derived for winged-circular and winged-elliptic cross sections (see 

appendix A). 

For a winged-circular cross section with th<* wing planform perpendicular to the cross- 
flow velocity (sketch (c)), we obtain, from equation (A 10), 


18 



V"o/ 

' 0 'Newt 


(2.27) 


For a winged-elliptic cross section with the semiir^or axis a and wing planform per- 
pendicular to the crossflow velocity (sketch (d)), we obtain, from equation (A 16), 




+ . 1 .- 1 

(bVfl*) a 


(2.28) 

For a winged-elliptic cross section with the semiminor axis b and wing planform 
perpendicular to the crossflow velocity (sketch (e)), we obtain, from equation (A21), 


/M =3^ 


a^lb^ 


[(aV**)- 1]®'^ 




n-' 


(2.29) 


»*’rom Newtonian impact theory, an expression for {CJC„ ) also has been derived 

" ”oNewt 

for winged-square cross sections with rounded comers (sketch (0). From equation (B1 5), 


hsH 


(^) 

°'Newt 


ir)k* ’ 


0<k<0.5 (2.30) 


where the equivalent diameter (from eq. (B12)) is 

(1 - (4 - ir)k* 




(2.31) 


Sketch (f) 

2.4.3 Values of (CjCf. ) and (CJC„ ) - From equations (2.23) through 

o SB o Newt 

(2.29), values of {CJC„ ) and (CJC„ ) were computed for circular and elliptic 

n ” ^oNewt 

cross sections alone and with wings. The results are plotted and compared in figures 7 
through 10. 

In figure 7, the variation of (C„/C„ ) with axis ratio alb is given for an elliptic cross 
section without wings. As previously noted in reference 9, values of (C„/C„ ) from 

n Rq 




slender-body theory are reasonably close to those from Newtonian theory for many ajb 
values of interest. 

In figure 8, the variation of with the ratio of wing semispans to body radius 

r is given for a winged-circular cross section. For s/r < 2, the values from both theories are 

reasonably close, but with further increase in s/r, the values of {CJC„ ) greatly exceed 

" "o SB 

those of (C /C ) 

”o Newt 

In figure 9, values of (C^/C^ ) are presented for a winged-elliptic cross section with the 
semim^or axis a perpendicular to the crossflow velocity V^. For the axis ratios of a/b = 2 
and 3, the figure gives the variation of {CJC„ ) with the ratio of semispans to semimajor 
axis a. As either alb or s/a increases, the disagreement between the results from the theories 
increases. 

In figure 10, values of (C„/C„^) are presented for a winged-elliptic cross section with 

the semiminor axis b perpendicular to the crossflow velocity I^. For axis ratios of a/b = 2 

and 3, the variation of {CJC„ ) with s/b is given. There is closer agreement between the 

‘ o 

values computed from the two theories for this cross-sectio ' 1 arrangement than for the 

arrangement where the semimajor axis and wing are perpendicular to V^. 

Table 3, taken from reference 9, shows the variation of (C„IC„ ) with comer 

" ^oNewt 

rounding k for square cross sections. The values of CJC„ computed from slender-body 

n Hq 

theory are reasonably close to those computed from Newtonian theory (eq. (2.30)). 


TABLE 3.- C„IC„ FROM NEWTONIAN THEORY FOR SQUARE CROSS 
SECTIONS WITH ROUNDED CORNERS 


k 

0 

0.05 

0.1 

0.2 

0.3 

0.4 

0.5 


1.33 

1.29 

1.25 

1.17 

1.11 

1.05 

1.00 


Corresponding values from slender-body theory vary from 1 . 1 9 at A: = 0 (no comer 
radius) to 1 .00 at k = 0.5 (completely circular cross section). 


2.S Relative Influence of Crossflow Terms 

It is interesting to examine briefly the relative influence of the potential and viscous 
crossflow terms in the equations for Cjy and C^. For demonstration, Jorgensen (ref. 9) 
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compared the magnitudes of the terms for an ogive-cylinder body of fineness ratio 1 1 
(ipj/d = 5) at a = 0° to 180° and A/oo = 2.9. The computed values (by eqs. (2.12) to (2.14)) 
are presented in figure 1 1. For this body, the viscous crossflow term contributes most of the 
normal force at high values of a and, of course, all of the normal force at a = 90°. Although 
the slender-body potential term contributes relatively little to at high a, it has a signif- 
icant influence on C^. 

The results in flgure 1 1 are indicative of those computed for most slender bodies 
(ref. 9). However, as shown in chapter S, the relative contributions of the crossflow terms 
can be modified considerably with the addition of thin lifting surfaces (wings and tails) to a 
body. 
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CHAPTER 3 

METHOD APPLIED TO BODIES OF REVOLUTION 

In recent years, there has been renewed interest in the basic aerodynamics of slender 
bodies of revolution because of emphasis on achieving more maneuverability from missiles. 
Some designs that use thrust-vector control systems are being considered for missile flight at 
angles of attack ranging from 0° to 180® (e.g., refs, li, 15, and 62). We will now assess the 
prediction method of chapter 2 by comparing predicted with measured longitudinal aero- 
dynamic coefficients for various bodies of revolution. 

3.1 Cone-Cylinder and Ogive-Cylinder Bodies at Af,, = 2.9 

We will compare predicted longitudinal aerodynamic coefficients with those measured 
by Jemell (ref. 63) for a series of three cone-cylinder and four ogive-cylinder bodies at 
a = -5“ to 180®. The free-stream Mach number is about 2.9 (Jemell quotes 2.86), and the 
Reynolds number based on body diameter is about 1 .25X 10* . The bodies with various nose 
and aftersection fineness ratios were tested in the NASA-Langley Unitary Plan Wind Tunnel. 

Figure 1 2 shows the seven bodies considered along with values of the geometric param- 
eten required to compute the aerodynamic characteristics. For the cone-cylinder bodies 
(numbers 1-3), all of the geometric parameters are easily computed, but for the bodies with 
tangent ogive noses (numbers 4-7) the required values of Apld^ , K* jd^ , xjd, and A^/d^ are 
not so easily obtained. Some convenient formulas (from ref. 9) for computing these param- 
eters for tangent ogives are given in appendix C. 

Equations (2.12) to (2.15), (2.18), and (2.19), along with the procedure outlined in 

chapter 2, have been used to compute the variation of Cyy, C^, and with a for the 

seven bodies considered. The values of Ca o and Ca . used in equations (2.18) 

^a=o "^a=i80 

and (2. 19) were computed in reference 9 with the assumption of turbulent boundary-layer 
flow. (Jemell (ref. 63) states that “boundary-layer transition was effected” by artificial trips 
throughout the a range.) The values include the contributions of fore pressure, base pres- 
sure, and the turbulent skin friction. The reader who wishes to make similar calculations is 
referred to reference 9. No attempt was made to estimate effects of wind-tunnel support 
interference. 
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In figures 13 to 16, computed values of Cj^, , and as a function of ot are 
compared with the experimental results for the seven bodies. Generally, there is close 
agreement of the computed with the measured results, especially in the variation of and 
Cffj with a. As expected, because of the approximation formulas used to predict with a, 
the poorest agreement is between the predicted and measured values of C^. 

It is satisfying to find that effects of afterbody fineness ratio, nose fineness ratio, and 
nose shape on and are predicted so well. Generally the magnitudes of and 
increase with an increase in fineness ratio, just as the computed results predict. Figure 13 
shows the effect of afterbody fineness ratio for cone-cylinder bodies, all with fineness- 
ratio-3 conical noses. Likewise, figure 14 shows the effect of afterbody fineness ratio for 
ogive-cylinder bodies, all with fineness-ratio-5 ogival noses. Figure 15 shows the effect of 
nose fineness ratio for ogive-cylinder bodies, and figure 16 shows the minor effect of change 
in nose shape from conical to ogival for a given nose fineness ratio of 3. 

In figure 17, computed positions of aerodynamic normal-force center (symbols) are 
compared with measured positions (lines from ref. 63) for the seven bodies at 2.9. 
The positions, are measured from the nose tip of each body in terms of the body 

length. As for the and results, the agreement of the computed with the measured 
values is reasonably close, especially for a near 90°. Note that symbols are used to denote 
computed values because only lines are given in reference 63 to denote the measured results. 

To assess the analytical method, plots of versus a may be omitted if plots of both 
Cjy and versus a are included. Hence, for conciseness, plots of versus a are omitted 
in the remainder of this study. Because precise prediction of versus a is beyond the 
scope of this study, plots of versus a are also omitted. 

3.2 Ogive-Cylinder Bodies at Mo^ = 0-6 to 2.0 

As shown in the previous section, the analytical method predicts the aerodynamic 
coefficients reasonably well for various cone-cylinder and ogive-cylinder bodies at a super- 
sonic Mach number of 2.9. We will now assess the method of predicting the and x^^ 
results for four ogive-cylinder bodies at a = 0° to 60° throughout the Mach number range 
from 0.6 to 2.0 (Moe - 0.6, 0.9, 1 .2, and 2.0). 
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Sketches of the four bodies considered are shown in figure 18. Bodies NjCi and NtCi 
were previously tested by Jorgensen and Nelson (ref. 16) in the NASA-Ames 6- by 6-Foot 
Wind Tunnel. Recently, these same bodies were retested in the same tunnel at = 0.6. 
They were also retested at Moo — 0-6. 0.9, 1.2, and 2.0 with cylindrical extensions attached 
to the aftersections. Cylinder C^, two diameters long, was attached to body fo form 
body N^CyCi, and cylinder C3, four diameters long, was attached to body N^Cx to fonn 
body N^C^C^. Thus, the bodies include nose fineness ratios of 2.5 and 3.5 (noses N-, and 
Nj), aftersection fineness ratios of 7, 9, and 11 (aftersections Ci, CiCj, and CiCs), and 
overall fineness ratios of 9.5 (/V7C1), 10.5 (iVjCi), 12.5 (NiCiCi), and 13.5 (N7C1C3). 

In figures 19-22, computed values of Cjy and (fi - as a function of a are com- 

pared with the measured results for the four bodies. The aerodynamic normal-force center, 
(6 - Xg^)ld, is measured forward on each body from its base in terms of its body diameter. 
For bodies N 2 C 1 and N7C1 (figs. 19 and 20, respectively), the results for all Mach numbers 
are compared on the same plots, and one can observe, at a given value of a, the variation of 
Cfj and {t-x^^)ld with Moo- For the longer bodies, N 2 C 1 C 2 (fig. 21) and N7C1C3 
(fig. 22), the results for each Mach number are compared on a separate plot. With these 
separate plots, there is less confusion in nparing predicted with measured values of 

(}enerally, the variation of with a is predicted closely for each body throughout the 
a and Moo ranges considered. At the lowest Mach number. Moo - 0*6, the Cj^ values previ- 
ously measured (ref. 16) for the shortest body, N-,Ci, are believed to be erroneous (fig. 20). 
The values from the retest agree closely with the predicted results. 

The aerodynamic normal-force centers are predicted best for the supersonic Mach 
numbers. CJenerally, for the subsonic Mach numbers, the predicted aerodynamic centers are 
more rearward on the boaies than the measured centers. For all Mach numbers and low 
values of a, the agreement of the predicted with the measured positions probably should be 
better than shown. The force and moment balance located inside each body was chosen to 
measure the large normal forces and pitching moments expected over the high a range. 
Thus, the accuracy of the experimental aerodynamic centers (determined from C^/C^) at 
the low values of oc is somewhat less than at the high values, and there is more scatter in the 
data. 
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In this study, we have concentrated primarily on predicting the static aerodynamic 
characteristics for bodies at high angles of attack. For bodies at low angles of attack, more 
theoretically exact prediction methods have been proposed. For example, for angles of 
attack less than about 20°, Perkins and Jorgensen (ref. 4) have shown that the agreement of 
theory with experiment at supersonic Mach numbers C5;n be improved by replacing slender- 
body potential theory (e.g., the first term in eq. (2.10)) with Van Dyke’s hybrid theory 
(ref. 64) or Tsien’s linearized theory (ref. 65). They also suggest the use of an empirical 
modification to the viscous crossflow theory (e.g., the second term in eq. (2.10)). In this 
modification, an experimentally determined correlation curve for the distribution of the 
crossflow drag coefficient along the body length is used when there is laminar crossflow. It 
is questionable whether the extra computation needed in these more detailed methods is 
warranted for most engineering studies. However, if an accurate loading distribution over a 
body length is desired, a more detailed method should be considered (ref. 4). 

3.3 Predicted Effect of Change in Crossflow Reynolds Number 
from Subcri'ical to Supercritical 

At present there is a general lack of aerodynamic data for which predicted results can 
be compared for bodies of revolution at supercritical crossflow Reynolds numbers and 
subcritical crossflow Mach numbers. Jorgensen (ref. 9), however, has computed results 
which demonstrate what can be expected for a slender body of revolution. The body chosen 
for study consists of a tangent ogive nose of fineness ratio 5 with a cylindrical aftersection 
of fineness ratio 6 (body 7 in fig. 1 2). 

For free-stream (and crossflow) Mach numbers less than about 0.4, the variation of 
and with a. has been computed for free-stream Reynolds numbers of 10* 10*, and 
10’’, and the curves are shown in figure 23. There is a significant effect of Reynolds number 
on both Cf^ and x^Ji throughout most of the a range. These curves, of course, reflect the 
strong influence of crossflow Reynolds number Re„ on crossflow drag coefficient Cj for 
two-dimensional circular cylinders (see fig. 2). 

As shown in figure 2, Cj for a circular cylinder drops considerably as Re„ increases 
from 10* (subcritical) to 10* (supercritical), and then there is a gradual rise as increases 
from 1 0* to 1 0’ (in the supercritical range). There is much more uncertainty in the magnitude 
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of at supercritical Re^ (such as 10* and 10’) than at subcritical Re„ (values less 
than about 2X10*), and the shading in figure 2 indicates the uncertainty because of scatter 
in known data. 

In figure 23, the shaded bands in the Cj^ and Xgji curves fot Re = 10* reflect the 
uncertainty in these curves resulting from the scatter in the Cj data shown in figure 2. It is 
clearly evident, however, that this uncertainty in the curves is relatively small compared 
with the large effect of change in Reynolds number. 

Figure 23(b) shows the ratio of Cj^ for the body at Re - 10* and 10’ to for the 
body at the subcritical Re of 10* . With this figure the effect of Re can be studied through- 
out the a range: for example, at a near 90°, the body at Re = 10* develops only about 
25 percent of the developed at Re = 10*, but at a < 10°, iOO percent of the C^is 
developed. Similar study of the movement of with change in Re can be made with the 
use of figures 23(c) and (d). 

Although no experimental data are available with which to compare these predicted 
effects for bodies of revolution, there are limited data for an early version of a noncircular 
space-shuttle body (ref. 66). Jorgensen (ref. 9) predicted some experimental trends (ref. 66) 
showing the decrease in with Of that results from an increase in crossflow Reynolds 
number from subcritical to supercritical. 
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METHOD APPLIED TO BODIES OF ELLIPTIC CROSS SECTION 
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We now compare predicted with measured normal-force and normal-force-center char- 
acteristics for three bodies of elliptic cross section at a = 0 ° to 60° and JUgo = 0.6 to 2 . 0 . 
First, we introduce the bodies considered, review the experimental test conditions, and 
present the equations used to compute and for each body. 

4. 1 Bodies Studied and Tests at = 0.6 to 2.0 

Figure 24(a) shows the three bodies of elliptic cross section considered here. Planform 
views of the bodies as they were oriented (in five different configurations) for the tests are 
shown in figure 24(b). All of the bodies were tested previously and the results presented in a 
data report (ref. 18). 

The basic circular body B, consists of a circular-arc tangent-ogive nose of fineness 
ratio 3 followed by a cylindrical aftersection of fineness ratio 7. Bodies and Bj have 
elliptic cross sections, ,md they have the same length and axial distribution of cross-sectional 
area as Bj . Hence the fineness ratio of ifd == 10 for B, is also the equivalent fineness ratio 
for Bj and B 3 , and all bodies have equal volumes. For Bj , the ratio of the semimajor to the 
semiminor cross-sectional axis (a /ft = 2) is held constant along the body length. Bodies B, 
and Bj were investigated in 1958 (ref. 6 ) only for a = 0° to 20® and Af^o = 2 to 4. Body B 3 
(new to the investigation in ref. 18) consists of the same nose shape as Bj but has an 
afterbody section of variable alb over four body diameters in length and a constant alb — 2 
over the rear three body diameters (see fig. 24(a)). Photographs of B 3 in figure 25 enable 
one to establish a clearer mental image of this more complex body. 

Six-component aerodynamic force and moment coefficients were measured in the 
study of reference 1 8 for these bodies in the Ames 6 - by 6 -Foot Wind Tunnel. However, in 
this study, only the variation of and (i - x^^)ld with a is considered. 

All bodies were tested at = 0.6, 0.9, 1 .2, 1 .5, and 2.0 and a = 0® to 58®. For the 
data used here, the Reynolds numbers, based on base diameter, are o.5X10* at = 0.6 
and 0.9 and 3.8X10* at Afoo= 1-2, 1.5, and 2.0. The elliptic bodies were tested at 0 = 0® 
(flattest side of nose pitching against the flow) and 90® (see fig. 24(b)). 
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4.2 Equations Used to Compute and for Each Body 

For the general case of a slender body in which the cross-sectional shape varies along 
the length, the normal-force and pitching-moment coefflcients can be computed from equa- 
tions (2.20) and (2.21): 

sin 2a cos 




and 


4 , 4 \ "of Newt 


C = 


sin 2a 




(4.1) 


27?C^^sin^a ^«/C„\ 


riXm - Jf)d>f 


(4.2) 


The axial distance from the body base to the normal-force center (see sketch (g)) is 
then given in terms of the body diameter by 


” X/iri ~ 


'ac’ _ m 

r* 

-N 


(4.3) 



28 



} 

1 




As previously discussed, (C„/C„ ) in equations (4.1) and (4.2) is the ratio of the 

"o SB 

local normal-force coefficient for the noncircular cross section to that for the equivalent 


(same area) circular cross section as determined from slender-body theory. The similar ratio 

(C„IC„ ) is determined from Newtonian impact theory. From equation (2.26), 

Newt 



( 4 .^^ 


where 0 is the angle of roll about the body longitudinal axis, being 0® with th.' semimajor 
axis a perpendicular to the crossflow velocity and 90® with the semiminor axis b perpen- 
dicular to the crossflow velocity (see sketches (h) and (i)). 



(h) (i) 

Sketches (h) and (i) 


From Newtonian theory, for the semim^or axis a perpendicular to the crossflow 
velocity 1^, equation (2.28) reduces to 




(4.5) 


For the semiminor axis b perpendicular to the crossflow velocity equation (2.29) 
reduces to 




(4.6) 


For the bodies of this investigation, equations (4.1) and (4.2) can be simplified by 
int nation. The simplified equations, which were used to compute C)y and values for 
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each body, are presented next. In all the and equations, empirical input values of n 
and Cj from chapter 2 are used. 

4.2.1 Equations for bodies with constant a/b cross sections (bodies Bt and B^)- Both 
bodies B\ and have cross sections of constant alb over the body length. For this 
condition, equations (4.1) and (4.2) reduce tc 


and 


where 



(4.7) 


(4.8) 


. (fn) 


4.2.2 Equations for body with variable a/b cross sections (body B^)- Body 3 has a 
midsection length of variable a/b, but the equivalent body of revolution isB, for which 
(L4/dx » 0 rearward of the nose-cylinder junction. Thus, equations (4.1) and (4.2) reduce to 



and 
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and Ah -Au> Also, for the bodies studied, 
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(4.10) 


Some potential or inviscid normal force can be expected intuitively to be carried over 
past the nose shoulder of body 5$, even though the approximate slender-body term gives 
zero normal force over the aftersection. For body at ^ = 90" (fig. 24(b)), the span 
increases with length past the nose shoulder, and it might seem reasonable to compute the 
slender-body terms for and based on the maximum span (or base) g^ treiry. F or 
this estimation. 
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- x)dx 
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4.3 Comparison of Con puted With Measured Normal-Force and 
Normal-Force-Center Characteristics 

In figures 26 through 28, computed valuer of and (C - x^^)!d versus a are com- 
pared with measured values (from ref 1 8) for the bodies studied. The comparisons are made 
for the bodies at a = 0° to 60° and free-stream Mach numbers of 0.6, 0.9, 1 .2, 1 .5, and 2.0. 

4.3.1 Bodies with constant ajb cross sections (bodies Bi and B 2 )— Generally, there is 
reasonably good agreement in figure 26 of the computed with the measured results for 
bodies Bi and B 2 (bodies of constant ajb along the length). The agreement is, however, 
better at the supersonic Mach numbers than at the subsonic. These comparisons, along with 
previous successful comparisons in chapter 3, tend to validate the prediction method as a 
useful tool in body aerodynamic studies, at least for bodies with circular and elliptic cross 
sections of constant ajb, 

4.3.2 Body with variable ajb cross sections (body B^)- As for bodies andi? 2 > the 

predicted characteristics for body B^ (with variable a/b) agree reasonably well with the 
measured results (fig. 27). However, the prediction of generally is not as close as for bodies 
Bi and £ 2 * The prediction of (C - ^ac^!^ ^ rem^^rkably close, especially at super- 

sonic Mach numbers. This close prediction of (C - x^^)jd, however, might be somewhat 
fortuitous since ihe prediction of is rot nearly as close. 

The prediction of with a is least accurate for body oriented at 0 = 90“. Gen- 
erally, for this case, is underpredicted at the lower values of a throughout the range. 
At the subsonic values of is underpredicted over the entire a range studied. 

Note that these predictions were made with equations (4.9) and (4.10), which do not 
account for any potential or inviscid lift rearward of the nose section. Intuitively and 
analytically (ref. 64), some inviscid normal force can be expected to be generated rearward 
of the nose. This is especially true for body 5, at 0 = 90° because the span increases with 
length past the nose, and it might seem reasonable to estimate the slender-body 
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potential-flow terms for and Cjjj based on the maximum span (see eqs. (4.11) and 
(4.12)). When the maximum span estimate is used, the prediction of with a improves 
somewhat (see fig. 28). However, the prediction of (£ - x^^)ld with oc deteriorates since the 
predicted values of (fi - x^^)ld move well forward of the measured values (fig. 28). 
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j CHAPTER 5 

I METHOD APPLIED TO BODY-WING AND BODY-WING-TAIL CONFIGURATIONS 

1 

? We now compare predicted with measured normal-force and normal-force-center char- 

acteristics for three body-wing configurations and three body-wing-tail configurations at 
a = 0° to 60° and Mgo = 0.6 to 2.0. First, we introduce the configuration models tested, 
review the experimental test conditions, and present the methodology used to compute % 
and for pach model. 

5.1 Configurations Studied and Tests = 0.6 to 2.0 

In references 20 and 21, many body-wing and body-wing-tail model combinations were 
tested. Planform views of the model components tested in the many model combinations 
are shown in figure 29. These components include a basic circular body Rj, an elliptic body 
B 2 with alb = 2, five flat-plate wings (R', to IVs), and a combination horizontal and vertical 
tail T. The bodies Bi and Rj are two of the bodies considered in chapter 4 (see fig. 24(a)). 
All the wings were designed to have the same planform area (16 d*) if the wings extended 
into the body to the axial centerline. Based on the phantom wing chord at the body 
centerline, the taper ratios for wings R', , W^, and were 0, 0.25, and 0.50, respectively 
(fig. 29(a)). They were also 0.25 for W4 and R'j (fig. 29(b)). Wings R', , W^, and R', 
(fig. 29(a)) had an aspect ratio of about 4; wings W^^, W^, and (fig. 29(b)) had aspect 
ratios of about 5, 4, and 3, respectively. Pertinent planform dimensions of the exposed parts 
of the wings are given in tables in figure 29. The tail dimensions are given in figure 29(a). 

Many combinations of these bodies, wings, and tail were tested (refs. 20 and 2 1 ) in the 
Ames 6- by 6-Foot Wind Tunnel at o: = 0° to 58° and — 0.6 to 2.0. The Reynolds 
number, based on body diameter d, was about 4X10®. 

Results from these tests (refs. 20 and 21 ) showed that changing wing taper ratio from 0 
to 0.5 changed the aerodynamic characteristics very little. Also, changing wing aspect ratio 
from 3 to 5 changed the aerodynamic characteristics very little. Thus, to achieve more 
conciseness in the present investigation, only three body-wing configurations and three 
body-wing-tail configurations were considered. 

; Planform views of the configurations studied in the present investigation are shown in 

|> 

f figure 30: the basic circular body Bi with wings of aspect ratio 3 and 4 (R^s and R'j, 
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respectively) and the elliptic body with the wing of aspect ratio 3 («'* ). These configura- 
tions (BiWs, B 1 W 2 , and B^W^) were also investigated with the tail T attached (configur- 
ations BiWsT, H'j T, and Bj M's 7^- 


5.2 Methodology Used to Compute and 
From equations (2.20) and (2.21), we obtain 
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and 
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r(x^ - .v)dx 


(5.1) 


(5.2) 


The first terms represent the Cjy and values from slender-body potential theory; 
the second terms represent the values from the viscous crossflow method modified by 
Newtonian theory. The first terms are not applicable, as v/ritten, for body-wing sections 
where the body (L4/dx values are zero or negative. Also, for body-wing and body-wing-tail 
configurations, more comprehensive methods from potential theory a.^ available. 

For this study, the first-term (potential) contributions to and were computed 
from the linear method presented in NACA Report 1307 (ref. 38). This method, referred to 
as the N-K-P method (for Nielsen, Kaattari, and Pitts), is restricted to bodies of circular 
cross section with wings and tails that do not have swept-forward leading edges or swept- 
back trailing edges. It is further restricted to small angles of attack and small angles of wing 
and tail incidence in which the forces are linear with angle. To obtain the wing-body 
interference, certain factors are defined that are the ratios of the lift on the components in 
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combination to the lift on the wing alone. These ratios are obtained primarily by slender- 
body theory, but the wing lift is obtained by linear potential theory. Wing-tail interference 
is treated by assuming one completely rolled-up vortex per wing panel and evaluating the 
tail load by strip theory. 

To combine the N-K-P method with the crossflow method, the N-K-P potential terms 
must be multiplied by a correction factor (sin 2a)/ 2a to produce a more correct type of 
nonlinear behavior to these terms and to eliminate the potential contribution as a 90° . 
With this modification. 


r* — //-* A j. 


J {%) ' 


rdx (5.3) 


and 
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m ff^'N-K-P 2a 


2r,C^^sin*a M ^c„\ 

J [%/ 


r(x^ - x)dx (5.4) 
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Since the N-K-P method is restricted to bodies of circular cross section with wings and 
tails, a further assumption must be made to estimate potential theory values of and 
for noncircular bodies with wings and tails. The local widths of the noncircular body in 


planform are replaced by the local diameters of a circular body, thus keeping the overall 


wing and tail spans constant (see sketch (j)). The crossflow method, of course, requires no 
such assumption for noncircular bodies. 



Assumed 
circular body 
for N-K-P 
method 


Sketch (j) 
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For the way in which the crossflow method is formulated in chapter 2, values of 

crossflow drag coefficient Cj for a circular cylinder are used. As shown in figure 1 , there is 

“w 

considerable variation of with crossflow Mach number for values of A/,, from about 
0.4 to 3. Also for less inai about 0.5, there is considerable change in as the 
crossflow Reynolds number Re^ exceeds the critical value of about 2X10* (figs. 2 and 3). 
These variations of with and may be expected for near-circular bodies, but 
surely not for very flat bodies or winged bodies. From the data available, values of Cj for 
flat bodies and plates do not appear to change nearly as appreciably with A/„ and Re^^ over 
the ranges shown in figures 1 and 2 (see, e.g., refs. 41, 66, and 67). Thus, for the body-wing 
configurations, it is likely that a constant value of Q will give closer agreement of theory 
with experiment, especially at a where is near or in the transonic regime. 

For flat-faced, two-dimensional configurations, reasonable values of crossflow drag 
coefficient can be computed from Newtonian or modified Newtonian theory (see, e.g., 
table 1 and ref. 66). For circular cylinders at low subsonic and hypersonic A/„, values of Cj 
computed from Newtonian or modified Newtonian theory also agree reasonably well with 
experiment (see fig. 1). In this study, modified Newtonian theory is used to compute the 
circular-cylinder Cj value that is substituted into equations (5.3) and (5.4). 

From modified Newtonian theory. 



- C 

3 ^Pstag 


1.2 forC, 


Pstag 


= 1.8 


For > about 4, Cp^^^^ss 1.8 from perfect-gas relations (e.g., ref. 68). In this study for 
wi '.g-body and wing-body-tail configurations, it is assumed that Cj =1.2 for all values of 
(and hence, Moo). 

The axial distance from the body base to the normal-force center is given (in body 
diameters) by 


^ " ^ac _ 1 ~ 

d C,v 


In this study, we present normal-force centers instead of C^„ values. 
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The entire method has been computer programmed for Ames Research Center by 
Nielsen Engineering & Research, Inc., Mountain View, Calif. The program is similar to those 
given in reference 69, and it combines essential parts of the N-K-P and crossflow programs in 
reference 69. It is written in Fortran IV for the CDC 6600 or 7600 machines. No tapes, 
drums, or disks other than the standard input/output units are required. The running time 
for a typical case is less than 1 sec. 

5.3 Comparison of Computed With Measured Normal-Force and 
Normal-Force-Center Characteristics 

In figures 3 1 through 36, computed values of and (fi - x^^)ld versus a are com- 
pared with measured values (refs. 20 and 2 1 ) for some of the body-wing and body-wing-tail 
configurations studied. As for the bodies alone, the normal-force coefficients are based on 
body cross-sectional area at the body base. Of course they can be easily converted to the 
more often used exposed wing planform area by dividing by A^lA^, where Ay^lA^ = 1 5.92 
for the bodies with and 16.49 for the bodies with 1^2 • The comparisons are made for the 
configurations at a = 0° to about 60°. For configuration IV 5 , - 0.6, 0.9, 1 .5, and 2.0 
(fig. 31). For the other configurations (fiiH'a, and B^WiT), 

Afoo = 0.6 and 2 . 0 . 

5.3.1 Body-wing configurations- Let us first look at the results in figure 31 for the 
circular body R, with the aspect ratio 3 wing Ws- At = 0.6 (fig. 31(a)), the measured 

values are predicted closely by the modified N-K-P potential method only for a up to 
about 10° or 15°. Then this potential method overpredicts for most of the higher a 
range considered. With even the potential method overpredicting Cyy, the combination of 
the potential and crossflow methods greatly overpredicts C^. It is to be expected, however, 
that at a near 90° the combination method should predict reasonably well, since Cyy is 
given entirely by the crossflow method (modified Newtonian impact theory) at a = 90°. 
Despite the difficulty in predicting over most of the high a range, the positions of the 
normal-force center are predicted closely by the combination method. 

With increase in Mach number to Afo© = 0.9, the comparisons of computed with mea- 
sured results are not significantly improved (fig. 31(b)). In fact, the normal-force centers are 
not predicted as closely as at A/^ = 0.6. Not until the free-stream Mach numbers become 
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supersonic do the modified N-K-P and combination methods give reasonable estimates of 
Cf^j. At Mao — I the measured results are only underpredicted a small amount by the 
modifi. d N-K-P method for a up to about 30° or 40° (fig. 31(c)). At A/oo= 2.0. the Ci^ 
results are predicted best by the combination method throughout the a range studied, 
although the experimental results are still overpredicted by this method (fig. 31(d)). It 
appears that the combination method should improve in its ability to predict as M^o 
increases throughout the supersonic range. 

With increase in wing aspect ratio from 3 to 4 (configuration to Bt M^j), the 

comparisons generally are not changed significantly. As shown in figure 32(a), for A/qo = 0.6 
the modified N-K-P and combination methods still overpredict Cyy significantly at high a. In 
fact, the overprediction for Bi W 2 is greater than for IV 5 (compare figs. 3 1(a) and 32(a)). 
As for the lower-aspect-ratio configuration, the results are predicted best by the com- 
bination method as the Mach number is increased to Moo ~ 2-0 (fig. 32(b)). 

When the body cross section is changed from circular to elliptical with the same aspect 
ratio 3 wing (configuration BxW^ to ^ 2 ^ 5 )’ th® results for Af<jo = 0.6 are predicted 
closely (possibly fortuitously) by the modified N-K-P method (fig. 33(a)). As before, when 
the Mach number is increased to Moo ~ 2 0, the combination method still gives the best 
prediction throughout the a range (fig. 33(b)). 

The break or significant deviation from linearity in the curve with increasing angle 
of attack at the subsonic Mach numbers makes it extremely difficult to formulate a rational 
method for predicting throughout the high a range. This break is attributed to flow 
separation over the wing upper surface. For wings of generally lower aspect ratio and higher 
leading edge sweep than those studied here, the break has been attributed to vortex break- 
down near the wing trailing edge of vortices shed from the leading edge. Vortex breakdown 
or bursting has been studied rather extensively by Wentz and Kohiman (ref. 70) for thin 
delta and nmdified delta wings with sweep angles from 45° to 85° at lew speed. They have 
observed that as a increases the position of vortex bursting of the trailing vortices moves 
upstream toward the trailing edge and crosses the trailing edge at a specific a. Above this a, 
a loss of lift occurs on the wing due to vortex bursting, and the effect becomes progressively 
larger as a increases. Mendenhall and Nielsen (ref. 71) have more recently collected data 
from several investigators for the a value at which vortex bursting occurs at the trailing edge 
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of delta wings tested at low speeds. They were unable to correlate the data and suggested 
that the factors which control vortex bursting were not reproduced or controlled between 
the various sets of test data. The wings used in the present investigation were generally 
swept less than those studied by Wentz and Kohlman (ref. 70) and Mendenhall and Nielsen 
(ref. 71). However, a 45'’ delta wing similar to that of Wi (fig. 29(a)) was investigated by 
^^entz and Kohlman (ref. 70). For this wing they failed to observe vortices, but they, of 
course, measured a loss in with increase in a over a particular value (near 20°). Despite 
the interesting research thus far, it seems that further research into the factors that control 
vortex bursting and flow separation is needed. The versus ct data presented here indicate 
that this is especially desirable for subsonic Mach numbers. Some initial exploratory flow- 
field pictures from recent vapor-screen and oil-flow tests are presented in chapter 6. 

5.3. 2 Body-wing-tail configurations- Now consider tne results for the same body-wing 
combinations but with the tail added (fips. 34-36). Generally, the comparisons of computed 
with measured results indicate similar trends as for the body-wing configurations. At 
= 0.6 and for the circular body fii , the results over most of the investigated a range 
are underpredicted, even by the modified N-K-P method (figs. 34(a) and 35(a)). With change 
to the elliptic body, however, the results (for T) are closely predicted up to a near 

50° (fig. 36(a)). At = 2.0, the results are predicted best by the combination method 

(modified N-K-P plus crossflow), although as for the body-wing configurations the measiired_, 
results are still overpredicted (fig. 36(b)). 

With the addition of the tail, there is generally more movement of the aerodynamic 
center (C - x^Jd) with a at = 0.6 (compare, e.g., results for Bi and B^ W^T in 
figs. 31(a) and 34(a)). This movement. v.i'ich takes place at a from about 10° to 40°, is only 
partially predicted by the combination method. It may be attributed to forebody and wing 
wake flow over the tail. Further investigation of the forebody and wing wake flow over the 
tail appears to be desirable. 

From the comparisons presented, it seems obvious that the methodology presented 
here represents only an initial step into the complex problem of predicting the aerodynamic 
characteristics of body-wing and body-wing-tail configurations to very high angles of attack. 
The reader interested in this field may wish to study several other initial approaches such as 
those of Mendenhall and Nielsen (ref. 7 1 ) and Axelson (re* 72). One should also include the 
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Polhamus suction analogy for wings (ref. 73) and some of its various adaptations and exten- 
sions (e.g., refs. 71, 74, and 75). Much additional research is necessary in the high a field, 
and this research initially should include visual observations of the flow fields. In chapter 6, 
we will show some of the photographs obtained from an exploratory visual study of the 
flow over bodies alone, bodies with a wing, and bodies with a wing and a tail. 


41 



1 


I 

i 


I 

I I 


CHAPTER 6 

VISUAL OBSERVATION OF FLOWS OVER MODELS 

It is a well-established fact that as the angle of attack (or incidence) of a model is 
increased above about 5“ or 10“ some flow separates from the upper surfp.ce, and vortices 
and/or regions of flow separation are formed. The basic vortex phenomena are easily 
described by referring to a body of revolution (see sketch (k)). At low to moderate angles of 
attack, two vortices are shed from the pointed body, and the vortex formation is sym- 
metrical (e.g., refs. 2, 5, 6, 13, 25, and 27). With further increase in angle of attack, the 
vortex formation may become somewhat asymmetrical. At some higher angle of attack, the 
feeding vortex sheets tear, and three or more vortices may appear (depending on the model 
geometry and frce-stream flow conditions as well as angle of attack). Both symmetric and 
asymmetric vortex formations are illustrated in sketch (k). Although generally illuminating 


Symmetric vortex formotion 




Asymmetric vortex formotion ond b.eokoway 


Sketch (k) 
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in concept, the illustrations in sketch (k) are somewnat simplistic, and further observations 
of vortex formations over many bodies alone and with wings and tails at various free-stream 
flow conditions are desirable. 

For observing the vortex formations over many models, the “vapor-screen technique” 
(e.g., refs, 2, 6, 76, and 77) is very useful and much less time-consuming than the more 
detailed pressure-probe-survey techniques (e.g., refs, 5 and 25). In this study, we have used 
the vapor-screen technique to observe the vortex formations over many of the same model 
configurations considered in chapters 3 through 5. In addition, we have added several cgival 
noses of different fineness ratio that attach to the body and body-win{ -tail combinations 
previously introduced. To observe surface flow and especially flow separation positions, the 
“oil-flow technique” (e.g., refs. 6, 17, and 76) also has been used for many of the same 
configurations. 

In this chapter, photographs are presented that show the flows over the various bodies 
alone and in combination with a wing and a wing plus tail at a = 10° i.o 50°. The free- 
stream Mach numbers are 0.6, 0.9, and 2.0, and the Reynolds number is about 4.3X10* 
based on body diameter. Before presenting the photographs, we will specify which models 
are considered and briefly review the vapor-screen and oil-fiow techniques used for the tests 
in the Ames 6- by 6-Foot Wind 1 unnel. 

6.1 Models Considered 

The models considered, as in the previous chapters, are identified according to body B, 
nose N, cylinder C, cylinder strake 5, wing W, and tail T. For the bodies of elliptic cross 
section (fig. 24), the roll orientation 0 is also specified. The configurations investigated and 
the figures in this report where the configuration dimensions are given are listed as follows; 


Configuration 

Figure 

= iV.C, 

24,37 

NiC, 

37 

N,C, 

37 


37 

NyCiS 

37 
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Configuration 

Figure 

0 = 0 ° 

24 

0 = 90° 

24 

B3 0 = 0° 

24 

B3 0 = 90° 

24 

B^Ws = A/.C.IV, 

29 


29 

B^W^T = N^C^W3 

T 29 

A/jC, W3 

29,37 

.VjCih'jT’ 

29,37 


Man-view sketches of these configurations are shown in figur: 38. 

Figure 37, which has not been presented before, shows a series of circular-arc tangent 
ogive noses that attach to the circular cylinder Cj . Noses A^, , A'j , and Nj have fineness 
ratios of 3, 3.5, and 5, respectively. Nose A /4 is formed by rounding the tip of a fineness- 
ratio-3.5 ogive (such as A/j ) to give a resulting fineness ratio of 3 . 

The thin wings and tail shown in figure 29 also attach to the cylinder C, , so many 
combinations of nose, cylinder, wing, and tail are ; >ssible. For this report ^nly some 
representative combinations are considered. 

6.2 Vapor-Screen and Oil-Flow Techniques 

6.2.1 Vapor-screen technique- In the vapor-screen technique, the wind tunnel is run 
with moist air. In fact, water is added to the airstream as needed. 

For supersonic Mach numbers, as the moist air expands through the supersonic nozzle 
into the test section it cools, and the moisture condenses to form a fog. This fog is illumi- 
nated by a sheet of bright light produced by high-intensity mercury-vapor lamps and pro- 
jected through the tunnel window(s) and across the stream. This sheet of light appears as a 
uniformly lighted screen of fog particles in the absence of a disturbance. However, with a 
model in the stream the uniform distribution of fog particles is disturbed, and the model 
disturbance affects the light scattered by the water particles. Wakes and vortices typically 
appear as dark “holes” in the screen. 
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For very high subsonic (A/oo = 0.9) and transonic Mach numbers, similar dark regions 
may appear at the vortex locations, but for lower subsonic flows = 0.6) light (conden- 
sation) regions often appear at the vortex locations. There are, in fact c '-^s where a light 
area will appear above a left wing panel and a dark area above the righi •>. and both areas 
will indicate separated regions and/or vortex flow. The pnysics may not always be simply 
explained, but the reader interested in more detail is referred to the treatise of McGregor 
(ref. 77). 

Figure 39(a) shows a schematic drawing of the vapor-screen apparatus for the Ames 6- 
by 6-Foot Wind Tunnel. Two light-source boxes and various camera locations are indicated. 
In this study, two light boxes were used for observations at = 0.6, but only one was 
required for necessary illumination at - 0.9 and 2.0. Each box contains six 900-W 
mercury -vapor lamps (BH-6) and mirrors to reflect the light through collimating slits. The 
light boxes (shown in fig. 39(b)) connect to a common shaft that passes over the tunnel, and 
the boxes can be pivoted in tandem and moved so that the light screen can cut the model 
longitudinal axis at the desired positions. A Honey well-Pentax spot light meter was focused 
on the vapor screen and moved with the light box shown in Hgure 39(b) to aid in maintain- 
ing consistent illumination for the photography. 

For the present investigation, a Hasselblad 70-mm still camera (node! 500 EL/7) was 
mounted on the sting support rearward of the model base (fig. 39(c)). The camera was 
enclosed in a protective housing. The gun camera shown mounted above the Hasselblad in 
figure 39(c) was not used for this study. All model.^ were painted black to minimize light 
reflection and to improve the quality of the photographs. The Hasselblad .amera had a 
50-mm f4 (wide-angle) lens, and TRI-X (400 ASA) film was used. 

6.2.2 Oil-flow technique- Flow patterns on the model surfaces at angles of attack 
were visualized through use of the oil-flow technique. In »his technique, the m» 'tHs were 
covered with a mixture of oil and titanium dioxide (TiO; ) and then run wet in the wind 
tunnel. The formula for the mixture was 5 teaspoons of SAE 30 oil and 5 teaspoons of 
TiO] , with about 3 drops of oleic acid added as an anticoagulant. To provide good contrast 
of the mixture with the models, all model' were fir A painted a flat black. Photographs of 
the oil-flow patterns on the models were taken during each run with Hasselblad cameras 
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focused through a side window and a special small window at the top of the test section of 
the Ames 6 - by 6 -Foot Wind Tunnel. 

6.3 Photographs Obtained from Vapor-Screen and Oil-Flow Techniques 

6.3.1 Photographs from vapor-screen technique Photographs taken with the vapor- 
screen technique for the models at = 0.6, 0.9, and 2.0 are shown in figures 40 to 54. 
Figure 40, presented as a prelude to the other figures, identifles some of the various items 
that appear in the other figures: the light sheet, light shadow, mcdel support, model base, 
vortex ••egions, vortex feeding sheets, and local shocks. The two photographs in figure 40 are 
for body N 3 C 1 at a « 30" and = 2.0. Light sheets are shown at the base of the nose 
(station 1 ) and at a more aft cylinder position (station 2 ). In the remaining vapor-screen 
photographs (figs. 41-54), a photograph of the flow field at the body base (station 3) is also 
included. The positions of the three flow-field stations are indicated on a model sketch in 
each figure. A model shadow similar to that shown in figure 40 usually does not appear on 
the photographs for - 0.6. As previously mentioned, at = 0.6, two light-source 
boxes were used in tandem, and light came from both sides of the wind tunnel. 

For a body of revolution (B, -NiC^) at Afoo' = 0.6 (fig. 41(a)), the simplistic flow 
model of the two rather tightly rolled-up vortices shown in sketch (k) is not evident, even at 
a < 30®. Rather there appear to be separation regions from both sides of the body that seem 
to almost coalesce along the lee side but do not roll up. In fact, t\/o narrow separation 
sheets very close together appear to trail back from the nose over the body length. With an 
increase in Of to about 40®, the separation regions from the nose appear more like the usual 
vortex regions. However, these vortex regions are very light in color in contrast to the dark 
regions at lower ot, and these vortex regions become asymmetric as they trail back over the 
body. Also, more than two regions develop. 

With increase in Mach number to = 0-9 and 2.0 (figs. 41(b) and (c)), the more 
traditional vortex formations appear, where there are two rather symmetric vortices (dark 
holes in the vapor sheet) shed from the body. Along both sides of the body, separation 
sheets “feed’ the vortices, and, of course, the vortices grow with movement from the nose 
to the base of the body. At = 2.0, local shocks from the vortex regions appear 
(fig. 41(c)). 
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From observation of figures 41 to 43, the effect of increasing nose fineness ratio 
from 3 (nose Ni) to 5 (nose 1V3) can be studied. Generally, with an increase in nose fineness 
ratio, the vortex formation becomes asymmetric at a lower a. Compare, for example, 
photographs for TV, C] ata»=30° (fig. 41) with those for N3C1 at o 30° (fig. 43). For all 
of these bodies, as for the remaining models to be discussed, vortex asymmetry is the worst 
at the lowest Mach number investigated (Mog = 0.6). Asymmetric flow separation and 
vortex asymmetry are accompanied by undesirable side forces (considered in chapter 7). 

When the tip of a fineness-ratio-3.5 nose (N 2 ) was rounded to make a fineness-ratio- 
3 nose (A^4 ), a strange vortex pattern developed at Mgg = 0.6 for the body N^Ci at a « 40° 
and 50°. This pattern (fig. 44(a)) consists of two very symmetric vortices from the nose 
located above two separation regions stacked one on top of the other (at station 2). The 
entire unusual pattern appears to be symmetric even up to a 50°. However, at the higher 
Mach numbers {Mgg = 0.9 and 2.0), the patterns are again similar to those forN2^i (com- 
pare, e.g., figs. 42 and 44). 

When strakes were attached to the side of the cylinder C, of body N^Ci, the vortex 
patterns became more symmetrical for a up to about 30° (compare figs. 43 and 45). How- 
ever, at about 40° and above, there was no effect. Apparently, the asymmetric pattern, 
originating with the fmeness-ratio-5 nose, could not be influenced by the strakes back on 
the cylinder. Note that the hand of the asymmetry of the vortices can be either left or right, 
and occasionally the pattern will switch even while it is being observed during a test run. 

When the body cross section was changed from circular (body 5i ) to elliptic (i ' B 2 
with alb = 2), the vortex patterns became more symmetric with the body oriented at 0° 

(flattest side toward the flow). This can be seen by comparison of figures 41 and 4 >w- 

ever, when the elliptic body B 2 was rolled to 0 = 90°, the vortex patterns became more 
asymmetric (see fig. 47). 

From tests of body B 3 , the body of elliptic cross section with variable a/b, the vortex 
symmetry was influenced mostly by the nose. As shown in figure 48, when Bj was oriented 
at 0 = 0° (flattest side of nose toward the flow), the vortex patterns were essentially 
symmetric at all test conditions. However, v/hen B 3 was rolled to 0 = 90° (thinnest side of 
nose toward the flow), the patterns became quite asymmetric at many test conditions (see 
fig. 49). 
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When wing it's of aspect ratio 3 was attached to body Ex , the body vortex growth over 
the wings was retarded, but extensive regions of How separation above the wing appeared 
(fig. 50). These regions were especially large and diffuse for = 0.6 and 0.9 (figs. 50(a) 
and (b)). For Mgg — 2.0, distinct vortices formed, originating at the forward wing-body 
juncture at the leading edge, on the wing upper surface, and at the wing tips. Thus, at the 
base of the body there were a pair of vortices from the nose, a pair from the wing-tip 
juncture, a pair from the wing leading edge, and a pair from the wing tips (see fig. 50(c)). 

In contrast with the flow at Mgg = 2.0, the flow at = 0.6 is very diffuse, and there is 
the possibility that the phenomenon of vortex bursting or breakdown (e.g.. refs. 70 and 7 1 ) has 
taken place at the higher angles of attack. This phenomenon (mentioned in chapter 5) might 
be indicated in some of the vapor-screen pictures by the light-colored diffuse vortex regions. 
Such regions can be observed at Mgg = 0.6 for both bodies alone and with wings, especially at 
the higher angles of attack. The phenomenon of vortex bursting might help explain the fact 
(demonstrated in previous chapters) that the characteristics for the bodies alone and with 
wings and tail are generally predicted best at the supersonic Mach numbers. The prediction 
methods, of course, do not account for this phenomenon. 

When wing of aspect ratio 4 was attached to body , similar flow patterns were 
obtained as with wing of aspect ratio 3. This can be confirmed by comparing the 
photographs in figures 50 and 51. Note, however, that the vortices shtd from the wing tips 
of Bi Ii ^2 3/qo = 2.0 lie outside the photograph frames. 

With the addition of the tail T to configuic^ion B^W^, there were no appreciable 
changes in the flow patterns (see photographs in figs. 5 1 and 52). 

When wing W 2 was attached to body N 3 C 1 , the nose-cylinder configuration that 
produced the greatest asymmetry at Mgg = 0.6 and 0.9, the vortex asymmetry from the 
fineness-ratio-5 ogival nose still persisted (compare figs. 43 and 53). This asymmetry also 
persisted when tail T was attached to the configuration (see fig. 54). 

6.3.2 Photographs from oil-flow technique- To obtain visual indications of the flow 
over the surface of the models, the oil-flow technique was used. Photographs from this 
technique support those from the vapor-screen technique in that flow-separation positions 
are quite clearly defined. As illustrated in a schematic of a body crossflow plane (sketch (1)), 
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oil ridges form near separation regions, and flow symmetry or asymmetry can be correlated 
with results from the vapor-screen technique. 



Sketch (1) 

Photographs taken with the oil-flow technique for some selected models = 0.6, 
0.9, and 2.0 are shown in figures 55 to 65. Both planform and side views are shown at the 
specified angles of attack. To avoid a tunnel support at the top of the test section, the 
planform views were taken with the camera somewhat off center, so the model planforms 
are not completely symmetrical. For convenience in comparing the oil-flow and vapor- 
screen photographs, the figures for each technique are indexed as follows; 


Configuration 

Oil-flow 

figures 

Vapor-screen 

figures 

B, = N,C, 

55 

41 


56 

42 

N3C1 

57 

43 


58 

44 

N 3 C 1 S 

59 

45 

Bi 0 = 0* 

60 

46 

Bj 0 =90° 

61 

47 
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Configuration 

Oil-flow 

figures 

Vapor-screen 

figures 

B 3 0 = 0° 

62 

48 

0 = 90° 

63 

49 

= NiCiWjT 

64 

52 

N^CyW^T 

65 

54 


The reader, of course, can make detailed observations of the comparable oil-flow and 
vapor-screen photographs for configurations of particular interest. Generally, when there is 
asymmetric flow separation from a model surface (indicated by the oil flow), there is also 
asymmetric arrangement of the vortices in the flow field (indicated by the vapor screen). 
Again, note that the hand of the asymmetry can be either left or right. Thus, because the 
oil-flow and vapor-screen photographs were not taken simultaneously, it is possible to 
observe a left-hand separation asymmetry and a right-hand vortex asymmetry or vice versa. 
The observer also should be aware that some local flow disturbances in the oil flow result 
from joints and wax-filled screw holes in the bodies. 
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CHAPTER 7 

EXPERIMENTAL SIDE FORCES ON MODELS AT /3 =.0° 

When models are pitched to high angles of attack, side forces can occur on the models 
even at zero sideslip angle (e.g., refs. 12-21, 78, and 79). These side forces generally occur at 
angles of attack between about 20° to 60° and in the subsonic-transonic Mach number 
range. They result from asymmetric flow separation and vortex flow over the leeward side 
of the models as shown in the photographs presented in chapter 6. 

As noted in a recent paper by Keener, Chapman, and Kruse (ref. 78), some aircraft 
have been lost due to uncontrolled flight at high angles of attack, and some of the loss in 
controllability might have originated from undesirable side forces and yawing moments 
attributed to flow separation and vortex asymmetry. Research on this phenomenon recently 
has increased considerably because the flight envelopes of modern aircraft and missiles aie 
being extended into the higher angle-of-attack range. 

We will now present and discuss some side-force data obtained recently (refs. 16, 
18-21) in the Ames 6- by 6-Foot Wind Tunnel for most of the model configurations 
considered in previous chapters (especially chapter 6) and shown in figure 38. For the bodies 
alone, we will discuss the effects of nose-fineness ratio, Mach number, nose-tip rounding, 
afterbody side strakes, and elliptic cross section. Then we will consider the effect of adding 
a wing and a wing plus tail to a body. 


7.1 Bodies Alone 

7.1.1 Effects of nose-fineness ratio and Mach number- In figure 66, the effect of 
nose-fineness iwUo 9.^/d on side-force coefficient Cy and center position (8 - x^^)ld for 
ogive-cylinder bodies of revolution is shown. Both Cy and (8 - x^^)jd (measured from the 
body base) are plotted as a function of angle of attack a for a up to 60°. Plots are presented 
for Mf)o = 0.6, 0.9, and 1.2. The models were tested at Mach numbers up to Mf„ - 2, but 
there were no side forces above about A/,jo= 1.2. The magnitudes of the side forces decrease 
with increasing Mach number. 

It is readily apparent from figure 66 that the largest side forces were obtained with the 
fineness-ratio-5 nose attached to the cylinder afterbody Cj (of fineness ratio 7). With 
decrease in nose-fineness ratio, the side forces decrease, and they almost disappear over most of 
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the a range for the body with the fineness-ratio-3 nose. For body N 3 C 1 (with the fineness* 
ratio-5 nose) at = 0.6, it has been shown that the maximum value of side force can 
become about 40 percent of the normal force (ref. 16). Keener et al. (ref. 78) have recently 
shown that the maximum value of Cy increases even more (to the order of or greater) as 
the Mach number is decreased from 0.6 to 0.25. There is no question that the side-force 
coefficients increase with decrease in Mach number and/or increase in nose-fineness ratio. 
Wardlaw and Morrison (ref. 79) also support these conclusions from their recent correlations 
of collected data. 

Note that the side-force centers seem to start well back on the cylinder; then as a 
increases above about 25°, they move forward onto the nose (see, e.g., fig. 66 ta)). At some 
much higher a (say about 55° to 60°) they then tend to move back onto the cylinder. Note 
also that for N^Ci , the body with the largest side forces, the maximum values of Cy are 
located well forward on the body. 

The side forces can be studied in conjunction with the vapor-screen photographs shown 
in chapter 6 . In figure 67, the data for N 3 C 1 at = 0.6 and 0.9 are plotted along with the 
vapor-screen photographs (from figs. 43(a) and (b)) taken at a = 10°, 20°, 30°. 40°, and 
50°. The vapor-screen photographs shown in figure-67 were taken at a crossplane station 
3.5 diameters forward of the model base. It is obvious that when the vortices become 
asymmetric (between a = 20° and 30°) the side forces develop. So long as the vortices 
remain steadily asymmetric, the side forces remain. Of course, as shown in chapter 6 , an 
asymmetric vortex system results from asymmetric flow separation from the body surface. 
There is obviously a need to analytically model the instability process that leads to asym- 
metric separation, asymmetric vortices, and undesirable side forces. 

In lieu of a theoretical method for computing the angle of attack for onset of side 
force, experimental procedures must be used. Keener et al. (ref. 78) have found that the 
onset angles can be roughly correlated with nose-fineness ratio (or semiapex angle) as 
families of curves of constant afterbody. They ha’-e found that the angle jf onset of side 
force for a given body is essentially invariant with Mach number, and, for a nose with no 
afterbody, the onset angle is given approximately by two times the semiapex angle. The 
general effect of increasing afterbody length is to decrease the angle of onset; that is, the 
longer the afterbody the smaller the angle at which a side force is first encountered. A 
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somewhat similar finding is reported by Wardlaw and Morrison (ref. 79). They also have 
produced (ref. 79) a rough correlation of data for the angle of attack at which the maximum 
side force is observed. It tends to decrease with increasing Mach number and body-fineness 
ratio. 

Keener et al. (ref. 78) also have made a rough correlation of data for the “upper-limit” 
angle of attack at which the static side force disappears and the wake flow becomes essen- 
tially oscillatory like a Karman vortex street from a two-dimensional cylinder. This upper- 
limit a also tends to decrease with increasing Mach number, varying from about a maximum 
of 80® at = 0.25 to a minimum of 50® at = 0 . 9 . 

7.1.2 Effect of nose-tip rounding- When the tip of a fineness-ratio-3.5 nose (Nj) was 
rounded to make a fineness-ratio-3 nose (A^ 4 ), the side forces at - 0.6 and 0.9 almost 
disappeared. This is shown in figure 68 where Cy and its center position are plotted against 
a for bodies Ni C\ , N2C1, and N4C1 . As discussed in chapter 6 , this nose rounding brought 
more symmetry to the flow field. The round-nosed body N4C1 , however, appears to be no 
better than the sharp-nosed body Ni Ci of the same fineness ratio (C^/</ = 3 ). 

7.1.3 Effect of afterbody side strakes- When strakes were attached to the sides of the 
cylinder Cj of body N^Ci (i^/d ~ 5), the side forces were not significantly changed. As 
shown in figure 69, the variation of Cy with a was changed somewhat, but the maximum 
values of Cy were about the same. 

7.1.4 Effect of elliptic cross section- The effect of elliptic cross section and roil angle 
on side-force coefficient and position is shown in figure 70. Results are compared for body 
Bx (a circular body of (C/d= 10 )) and body (the equivalent elliptic body of constant 
alb = 2) at Mfg, = 0.6 and 0.9. With the elliptic body oriented at 0 = 0® (flat side 
pitching against the free-stream flow), the side-force coefficients are very small and close to 
those for 5,. However, when B 2 is rolled to 0 = 90®, the side-force coefficients increase 
considerably. According to reference 18, Cy becomes more than twice at a = 50° for 
B 2 at 0 = 90®. As shown in chapter 6 , the separation and vortex patterns were very sym- 
metric for Bj at 0 = 0 °, but they became very asymmetric for Bj at 0 = 90®. 

From tests of body Bj, the body of elliptic cross section with variable alb, it was 
found that the wake flow-field asymmetry and side forces were influenced mostly by the 
nose. As shown in figure 7 1 , when B 3 was oriented at 0 = 0® (flattest side of nose pitching 
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against the free-stream flow), the side-force coefficients were very small and close to those 
for bx (the equivalent body of revolution). However, when was rolled to 0 = 90°, 
undesirable side forces developed, and, as shown in chapter 6, the vortex flow was quite 
asymmetric (fig. 49). 

As for the bodies of revolution, the maximum values of side-force coefficient decreased 
as the free-stream Mach number increased. At = 1-5, there were essentially no side 
forces measured on any of the elliptic bodies even for the worst roll orientation (ref. 18). 

7.2 Body-Wing and Body-Wing-Tail Configurations 

7.2.1 Effects of adding a wing and a wing plus tail to a body - In figure 72, the effects 

on side-force coefficient and side-force position of adding a wing and a wing plus tail to a 
body are shown. The body A^sCi, which has a fineness-ratio-5 nose, produced the largest 
side forces for the bodies of revolution tested. As shown in figure 72, these side forces and 
their positions remain about the same with the thin wing aspect ratio 4 attached or 

even with the tail T added. (Dimensions for the wing and tail are given in fig. 29 and for the 
body in fig. 37.) 

From these comparisons it can be concluded that the most important influence comes 
from the body nose. As mentioned in chapter 6, the vortex asymmetry that appeared from 
the fineness-ratio-5 nose, when tested with only the afterbody '’yiinder C, , still persisted 
when the wing and the wing plus tail were attached (see, e.g., figs. ' ’ and 53). 

7.2.2 Effects of wing-aspect ratio and taper ratio- For thin wings of about equal 
planform area (fig. 29) but with aspect ratios of 3, 4, and 5, there was essentially no effect 
of aspect ratio on the measured side forces of the body-wing models tested (ref. 20). 
Likewise, a change in taper ratio from 0 to about 0.5 resulted in no appreciable side-force 
effect (ref. 20). The results were also unchanged when the tail T (fig. 29) was attached to 
the wing-body models (ref. 21). 
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CHAPTER 8 

CONCLUDING REMARKS 

A review and an extension of an engineering-type method have been presented for 
computing the normal-force and pitching-moment coefficients for slender bodies of circular 
and noncircular cross section alone and with lifting surfaces. In this method, a semiempirical 
term representing viscous separation crossflow is added to a term representing potential- 
theory crossflow. In the generalized equations written for Cjy and C^, ratios are required of 
the local normal-force coefficient per unit length for the cross section of interest so that for 
the equivalent (same area) circular cross section. These ratios are given both from slender- 
body and Newtonian theories. Formulas and numerical values of thest ratios are included 
here for winged-elliptic and winged-square cross sections, the square cross sections having 
rounded corners if desired. 

In computing normal-force and pitching-moment coefficients for the bodies alone, 
slender-body theory was used for the term representing potential crossflow. In computa- 
tions for the bodies with thin wings and a tail, the linearized potential method of Nielsen, 
Kaattari, and Pitts was used, modified for high angles of attack. 

For many bodies of revolution, computed aerodynamic characteristics were found to 
agree well with measured results for investigated frec-stream Mach numbers from 0.6 to 2.9. 
The angles of attack ranged from about 0° to 180“ (or = 2.9 and from about 0“ to 60° 
for Mgg = 0.6 to 2.0. Agreement of predicted with measured results was best at supersonic 
Mach numbers. 

For several bodies of elliptic cross section, measured results were also predicted reason- 
ably well over the investigated Mach number range from 0.6 to 2.0 and at a = 0“ to 60°. As 
for the bodies of revolution, the predictions were best for supersonic Mach numbers. The 
predictions were better for a body in which the cross-sectional shape (alb) remained con- 
stant over the length than for a body in which it varied. Although the prediction technique 
probably can be improved with further research, it is felt that the predictions are suffi- 
ciently accurate for most preliminary design studies. 

For body-wing and body-wing-tail configurations with wings of aspect ratio 3 and 4, 
measured normal-force coefficients and normal-force centers v^’ere predicted reasonably well 
at the upper test Mach number of 2.0. However, with a decrease in Mach number to 


55 



t 

I 

I 


/ 


\ 


M^ = 0.6, the agreement for Tyy rapidly deteriorated, although the normal-force centers 
remained in good agreement. At the subsonic Mach numbers = 0.6 and 0.9) and angles 
of attack above 10° or 20°, the measured results were even overpredicted somewhat by just 
the potential-flow term of the combination method. 

From vapor-screen and oil-flow studies, it was observed that the How was completely 
separated over the wings at the subsonic Mach numbers for a > about 20°. lliere was some 
evidence of vortex bursting breakdown for some of the models (including bodies alone) 
at subsonic Mach numbers. At Mgo = 2, however, discrete vortices from the bodies and 
wings were observed at the base of the models. 

For many body, body-wing, and body-wir>j-tail configurations, vapor-screen and oil- 
flow photographs were obtained for 0-6, 0.9, and 2.0 and a - 10°, 20°, 30°, 40°. and 
50°. It has been observed that, when the separation and vortex patterns were asymmetric, 
undesirable side forces could be measured on the moueis even at zero sideslip angle. These 
side forces generally originated at a > 20°. 

For bodies alone, the side forces can be significantly affected by changes in Mach 
number, nose-fineness ratio, nose bluntness, and elliptic cross section. The side-force coeffi- 
cients decrease with increase in Mach number through the subsonic-transonic range and 
disappear with increasing Mach number into the supersonic range. 

From tests (at 0.6 <Afoo ^ 2.0) of tangent-ogive noses connected to a circular cylinder 
of fineness ratio 7, it has been found that the side-force coefficients increase from about 
zero for a fineness-ratio-3 nose to a maximum of about 40 percent of the normal-force 
coefficient for a fineness-ratio-5 nose. Other researchers have reported even more increases 
in Cy (of the order of or greater) as the Mach number is decreased from 0.6 to 0.25. 

Nose-tip rounding significantly decreased the side-force coefficients for a circular body 
with a sharp-nosed ogive of fineness ratio 3.5. However, the beneficial decrease was no 
greater than that obtained by merely using a sha*-p-nosed ogive nose of the same fineness 
ratio (fineness ratio 3) as that for the resulting blunted nose. 

When strakes were attached to the sides of the afterbody cylinder of an ogive-cylinder 
model with an undesirable fineness-ratio-5 nose, the side forces were not significantly 
changed. Thus, the influence of the nose was dramatized. 
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Gianges in cross section from circular to elliptic for an ogive-cylinder body (of 
fineness-ratio-3 nose and fineness-ratio-7 afterbody) produced some interesting effects on 
side-force coefficient. At all test Mach numbers (0.6 to 2.0) and angles of attack (0° to 60°), 
the side-force coefficients were generally small or negligible for the circular body and the 
generated elliptic bodies at ^ = 0° (flat side pitching against the stream crossflow). With the 
elliptic bodies at 0 = 90°, however, some values of side-force coefficient became as large as 
twice the values of normal-force coefficient at the same high angles of attack. From the 
standpoint of reducing undesirable side forces at high angles of attack, it was always found 
best to have the flattest side of the elliptic body nose pitching against the stream crossflow, 
even when ajb was not constant over the body length. 

Undesirable side forces measured for body-wing and body-wing-tail models were gen- 
erally about the same as those measured for the bodies alone. As for the bodies alone, the side 
forces developed at subsonic Mach numbers for a > about 20°. Also, as for the bodies alone, 
the side forces and yawing moments increased with increase in nose fineness ratio. Fineness 
rati<»s greater than 3 produced the largest side forces. No effects of wing-aspect ratio or 
taper ratio were observed. From these comparisons, it can be concluded that the most 
important influence comes from the body nose. 
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APPENDIX A 

DERIVATION OF (CJC,, ) FOR WINGED-CIRCULAR 

" "osewt 

AND WINGED-ELLIPTIC CROSS SECflONS 


Winged-Circular Cross Section with Wing Planform Perpendicular 
to Crossflow Velocity 


For blunt configurations of the type 
shown in sketch (m), it is assumed that each 
elemental particle of fluid strikes the 
configuration at velocity I^ and thereupon 
loses its normal component of momentum. 
Ihis leads to the well-known Newtonian 
expression for the pressure coefficient: 

Cp = 2sin^6 (Al) 



y 



Sketch (m) 


where 5 is the local angle a tangent to a forward-facing surface makes with the free-stream 
direction (sketch (m)). In Newtonian flow, the pressure coefficients over the rearward face 
are assumed to be zero. The total section crossflow drag joefficient 'based on body diam- 
eter </ = 2r) is then given by 



(A2) 


For a circular section as shown in sketch (m), 
.V* + = r* 


tan 6 


dv 



(A3) 

(A4) 
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sin* 6 = — 


I + tan* 6 


Equations (A 1 ) and (AS) are substituted into equation (A2) to give 


= 7/ 


C„ =2 sin* 6 = 2 
rwtng 


From integration of equation (A6), we obtain 




where, for the equivalent circular section, 

J (A9, 

Thus, the ratio of tlie normal-force coefficient per unit length for the winged-circular 
section to that for the equivalent circular section is given by 


iS) -f) 


Winged-Elliptic Cross Section with Semimajor Axis and Wing Planform 
herpendicular to Crossflow Velocity 

For a winged-elliptic cross section where the semimajor axis (a) and the wing planform 

are perpendicular to the crossflow velocity (see sketch (n)). the same procedure is 

followed In the derivation of ^CJCf. ) as introduced for the winged-circular cross 

o Newt 

section. Both the basic equations, (Al ) and (.\2>. are used again. 


,S9 



1 

L 


V 


n 



^2 


For an elliptic body as shown in sketch (n), 

(All) 




and 
sin^ 6 = 


2 f. _ tan^ 6 


- y'^ 


1 + tan^ 5 - y‘[ \ - (b^la^)\ 


(A13) 


oketch (n) 


With equation (A13) combined with equations (Al), (A2), and (A7), we obtain the 
section crossflow drag coefficient: 


c =i r ~ y"" 

^ dj - y‘^[\ - (b'^ la^)\ 


dy + j 




d.v 


(A 1 4) 


o a 

where d is the diameter o. the equivalent-circular body and is given by 

d = 2V^ 


(A15) 


Also, as for the winged-circular cross section, the crossflow drag coefficient for the equiva- 
lent circular cross section is given by = 4/3, 

Thus, by integrating equation (A14), substituting 2>/oF for d, and dividing the result 
by 4/3, we obtain 


(S) -ifi 

\ "o/.Vewr 


-b^ la^ 


[1 - 


log 


[f (' " l/^) 


-1- ^-5—;,- + — - i; 

1 - (b^la ) a 


(A 16) 
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since 


(t) i^) 

^ (^'Newt ^ ^'Newt 


Winged-Elliptic Cross Section witn Seminiinor Axis and Wing Planform 
Perpendicular to Crossflow Velocity 


For a winged-elliptic cross section 
where the semiminor axis b and the wing 
planform are perpendicular to the cross- 
flow velocity (sketch (o)), the same 
procedure is followed as for the previous 
configurations. 

For an elliptic body as shown in 
sketch (o), 


Vn 



= j 


a 


dx a* \y / a T 


and 


sin^ 6 - 


tan* 6 


i’^ 


1 -I- tan* 5 6* - >’*[1 - (a^lb^)] 


(A 17) 


(A18) 


(A19) 


With equation (A19) combined with equations (Al), (A2), and (A7), we obtain 

•b 


'd d 


r 6 * 

J b^-y^ 


**_-y^ H , + 1 

[r-(a*/ft*)l ■' d 


r. 


(A20) 


where d = 2y/aB is tb ' diameter of the equivalent circular cross section. As before, the cross- 
flow drag coefficient for the equivalent circular cross section is given by Cj ~ 4/3. 

Thus, by integrating equation (A20). substituting 2>/5F for d, and dividing the result 
by 4/3, we obtain 
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a^lb^ 

[(a^lb^)- 1 ] 


3/2 


tan ‘ 



since 


1 + ±- ]\ 

(a^lb^)- 1 b I 

(A21) 
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APPENDIX B 


DERIVATION OF iCjC„ ) FOR WINGED-SQUARE CROSS SECTIONS 
" "o Newt 

WITH ROUNDED CORNERS 


For winged-square cross sections with 
rounded comers (sketch (p)), the pressure coef- 
ficient over the front face is given by 


Cp = 2 sin^ 8 


where 6 is the local angle that a tangent to a 
forward-facing surface makes with the free- 
stream direction. In Newtonian flow, the pres- 
sure coefficients over the rearward face are 
assumed to be zero. The total section drag coef- 
ficient (based on width w) is then given by 


r=kw 


I 


Sketch (p) 


o 


2 

w J ^Pl 

o 

2 /•(W2)-r 


^Pbody w f Sv 


dy + - j C- dy + -^ I C„ iy (B2) 
Pflat vv I Pcorner w J Pwirtt, 


For a rounded corner, as shown in sketch (p). 


X* -t- y* = r’ 


,a„5 = $=- * =- 


sin* 8 = 


tan* 8 
1 -I- tan* 6 


- 1-4 



The pressure coefficient for the rounded corner, obtained by substituting equation (B5) into 
(Bl),is 


'corner 



(B6) 


The pressure coefficients for the front flat portion of the body and for the wing are 
given by 


Cp = 2 sin^ 6 = 2 

since 6 = 90° 

Thus, with equations (B6) and (B7) substituted into equation (B2), we write 


o = - 


4 /•(w/2)-r 


w 


/ 


dy - 
w 


o o 

From integration of equation (B8), we obtain 

4 


0 w/2 




(B7) 


(B8) 


(B9) 


Now, to obtain {CJC„ ) , it is necessary to find the diameter d of the equivalent 

o Newt 

circular cross section which has the same area as the cross section studied. The area of the 
equivalent circular cross section is, of course, 

A= (BIO) 

4 

and the area of the cross section studied is 

/4 = w*[l - ^*(4 - 7T)] (BID 


By equating equations (BIO) and (B1 1), we obtain the equivalent diameter 
With Cj from equation (B9) based on d instead of w. 


(B12) 


(B13) 
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b 

i 


where, for the equivalent circular section. 



Thus, from equations (B12), (B13), and (B14), we obtain 


(B14) 


( 1 ^) -[w) 

^ ^'Newt ' 


it)k} 


(B15) 


0 < it < 0.5 
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APPENDIX C 

FORMULAS TO COMPUTE GEOMETRIC PARAMETERS FOR TANGENT OGIVES 


To compute the aerodynamic characteristics of bodies of revolution having tangent 
ogive nose shapes, various geometric parameters must be obtained. To compute Cyy and 
the planform area Ap is required; to compute C^, it is also necessary to obtain volume V 
and distance from the nose vertex to the centroid of planform area. To compute skin- 
friction drag, the wetted surface area is needed. 

For an ogival nose of length and diameter d (see sketch (q)) the following useful 
formulas have been derived (lef. 9): 



Sketch (q) 

_ 6^ (2/3)j/?* - l/?2 - \R - (1/2)1 

T ~ d ~ Apfd^ 

and 

where R is the ratio of the ogival arc radius to base diameter d and 



(C3) 


(C4) 


(C5) 
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Figure 2.- Variation of crossflow drag coefficient with crossflow Reynolds number for circular cylinders 

at subcritical crossflow Mach numbers (M,, < 0.4). 







Figure 3 - Variation of crossflow drag coefficient with crossflow Reynolds number for 
circular cyUnders at supercritical Reynolds numbers and at crossflow Mach numbers from 
0.25 to 0.50 (from ref. 54). 



Circular cylinder length-to-diameter ratio 


Flat plate length - to -width ratio 

Figure 4.- Ratio of crossflow drag coefficient tor a finite-Iergth cylinder (or flat plate) to 
that for an infinite-length cylinder (or flat plate) (from ret. 56). 
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BODIES Rp y in 

N.C. N,Ci 


Mn = Mflo Sin a 

>ure 5.- Variation of rj with )W,j (crossflow Mach number) computed from experi 
mental data (ref. 1 6) for bodies of revolution at high a (45° < a < 60° ). 

o 7] COMPUTED FROM FIGURES I AND 5 
a n FROM FIGURE ^ FOR i/d = 10 
o 77 FROM FIGURE 4 FOR i/d = 12 


Mn'Moo Sin a 

Figure 6.- Variation of tj with obtained from experimental results 
in figures 1,4, and 5. 
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(a) Semimajor axis a perpendicular to crossflow velocity 
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(b) Semiminor axis b perpendicular to crossflow velocity 


Figure 7 - Ratio of local normal-force coefficient for an elliptic cross section 
to that for the equivalent circular cross section. 
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Figure 8.- Ratio of local normal-force coefficient for a winged-circular cross section 
to that for the equivalent circular cross section. 
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(a) a!h = 2 



(b) a!h = 3 

Figure 9.- Ratio of local normal-force coefficient for a wingeii-clliptic cross section t > that 
for the equivalent circular cross section: seminiajor axis a perpendicular to crossllow 
velocity 
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s/b 

(b) alb = I 


Figure 10.- Ratio of local normal-force coefficient for a winged-elliptic cross section to that 
♦"or the equivalent circular cross section: scmin inor axis h perpendicular to crossflow 
velocity 
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Slender -body potential theory 

Slender-body potential plus 

VISCOUS cross flovi/ theory 



Slender body 



Viscous crossflow 


0 20 40 60 80 100 120 140 160 180 

a ,deg 

Figure 11.- Computed components ot norinal-toree and pitching-nu)ment coellicients lor 
an ogivc-cylinder bod’^ at a = 0° to 180° ;.nd - Urom ret. 9). 








I igure I 3. Comparison ot computed with measured aerodynamic characteristics tor cylindrical 
bodies with conical noses of tineness ratio 3: = 2.9. 
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(b) Cylindrical bodies with ogival noses 
of fineness ratio 5. 


0 20 40 60 80 lOu 120 140 160 180 

a, deg 

(c) Cylindrical bodies with conical and 
ogival noses, overall fineness ratio 9. 


Figure 1 7.- Comparison of computed with measured aerodynamic force centers for cylindrical 
bodies with conical and ogival noses; 3/oa, = 2.9. 
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Figure 18. - Ogive-cylinder bodies for which the aerodynamic characteristics have 
been computed and measured for 0.6 < < 2.0. 
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Figure 21 . - Continued 






Figure 21. Continued 







Figure 21. - Concluded 
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planform and end vews of bodies 
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; (a) Modei dimensions with d = 6.6 cm (2.6 in.) 

Figure 2^.- Bodies ol elliptic cioss section for which the aerodynamic characteristics were measured 

in reference 18 and computed in this study. 
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.AN VIEW 



Figure 26.— Comparison of computed with measured aerodynamic characteristics for bodies 

B\ and B>i \ 0.6 ^ ^ 2.0. 






Figure 26.- Continued. 
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26.— Concluded. 



COMPUTED MEASURED 



Figure 27.- Comperison of computed with measured aerodynamic characteristics for bodies 

Bi andflj;0.6<M«,<2.0. 
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Figure 27.— Continued. 
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Figure 27.- Continued. 
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Figure 27.- Continued. 
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Figure 27,— Concluded. 



C0NF1GU»ATICN COMPUTED MEASURED 



Figure 28.- Comparison of computed with measured aerodynamic characteristics for body 

S. at ^ = 90";0.6 2.0. 






(b) A/«, = 0.9, Rc = 6.5X10*. 
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Figure 28.— Concluded. 
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(a) Body Bi with aspect-ratio 4 wings of various taper ratios 
and tail arrangement. 
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(b) Body Bi with wings of aspect ratio 3, 4, and 5 (W ^ , W 2 , and ). 

Figure 29.- Components for body-wing and body-wing-tail models tested in 

references 20 and 21. 
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(c) Body Bi with wings IVj , K'a , and M's . 
Figure 29.- Concluded. 




Figure 30.- Planform views of configurations for which the aerodynamic characteristics 
were measured in references 20 and 21 and computed in this study. 
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Figure 3 1 Comparison of computed with measured aerodynamic characteristics for 

body-wing model B^W^. 
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Figure 31.— Continued. 
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Figure 31.- Concluded. 




Figure 32.— Comparison of computed with measured aerodynamic characteristics for 

body-wing model B^W^. 



} 

I 



o 

Ul Q. 

± 


a 

Ul 


8 % g 


I 

»s 

8 

oc 

o 

U. I 


ii§ 

I 


I 

I 




Figure 32.- Concluded. 






Figuie 33.— Comparison of computed with measured aerodynamic characteristics for 

body>wing model Aj It's . 
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Figure 33.- Concluded. 






Figure 34.- Comparison of computed with measured aerodynamic characteristu/; for 

body -wing-tail model Bt M's 7*. 
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Figure 36.- Comparison of computed with measured aerodynamic characteristics for 

body-wing-tail model M's T. 
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Figure 36.- Concluded. 





Figure 37.- Additional tangent ogive noses and strake for modification of body Bi . 



Figure 38.- Plan-view sketches of configurations used for vapor-screen and 

oil-How studies. 
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(b) Picture of light*box installation. 
Figure 39.- Continued. 
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(c) Picture of model and camera installation. 
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Kigure 39.- Concluded. 
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MODEL N^Ci 


STATION I — 


STATION 2 - 
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Figure 40. - Identification of typical Itema shown in vapor-screen photographs taken for 

body at * 2.0. 
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Figunr 41 . Vapoi-scrccn photographs for B, = .V, C’l (body with fineness-ratio-3 nose) 
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Figure 41 . - Concluded. 
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Figun.* 42. - Vapor-scnrcn p’lotographs forNjCi (bcidy with fineness-ratio-3.5 nose). 
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Figure 42.- Continued. 
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Figure 43. - Vapor-screen photographs for A^ 3 C, (body with fineness-ratio-5 nose). 
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Figure 43.- Continued. 
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Figure 44.- -screen photographs for N^Cx (body with blunt nose of fineness ratio 3). 
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Figure 44.— Continued. 
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Figure 45. Vapor-screen photographs tor .VjCii' (body with fiiieness-ratio-5 nose and afterbody strake). 
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Figure 45. Concluded 
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Figure 46. Vapor-screen photographs for B 2 at ^ = 0" (body with constant ajh cross sections). 
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Figure 47. - Vapor-screen photographs for ^ - ^0° (body with constant a/h cross sections). 





Figure 47. ~ Concluded. 
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Figure 48.- Vapor-screen photographs for Zf, at ^ = 0° (body with variable a/Z> cross sections). 
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Figure 48.- Continued. 
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F'igurc 49. Vapor-screen photographs for ^3 at 0 = 90° (body with variable a/h cross sections). 
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Figure 49. Concluded 
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l-ipurc 5 1 . V'apt'r-scrccn photographs lor It', = .V,C, h’, (bixly with finciicss-ratio-3 nt se and aspcct-ratio-4 wtngi. 
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Figure 52.- Vapor-screen photographs for B 1 W 2 T — iV, C| T (body wi^h fineness-ratio-3 nose, aspect-ratio-4 wing, and tail). 
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Figure 52. Continued. 
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Figure Vapor-screen photographs for /V 3 CJ (^50dy with fineness-ratio-5 nose and aspect-ratio-4 wing). 


STATION I 



I 76 


Figure 53. Continued. 
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Figure 54. - Continued. 
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Figure 54. Concluded. 
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Figure 56.- Continued 








(a) Maa - 0 - 6 . 

Figure 57.- Oil-flow photographs for N 3 C, (body with fineness-ratio-5 nose) 





(b) Afoo = 0.9. 
Figure 57.- Continued 
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(a) M^^O.6. 

Figure 60.- Oil-flow photognpht for Sj at 4 * O'* (body with constant a/b croM sections). 
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(a) Af«, ■ 0.6. 

Figure 62.- Oil-flow photographs for at 0 = 0* (body with variable alb cross sections). 
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PAGE IS Figure 62.- Continued. 
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(c) ^00 = 2.0. 
Figure 62.- Concluded. 
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(a) A#oc * 0.6. 

Figure 64.- OiFflow photographs for BiW^T • C| T (body with fineness-ratio*3 nose, 

aspect*ratio4 wing, and tail). 
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(a) ,W<* = 0.6 - Concluded. 
Figure 64.- Continued. 
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(a) Moo = 0.6. 

Figure 65.- Oil-flow photographs for N 3 C 1 W 2 T (body with fineness-ratio-5 nose, 

aspect-ratio-4 wing, and tail). 
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(a) Moo - 0-6 - Concluded. 
Figure 65.- Continued. 


215 







f 


IDE 


K»T-mJN 
SAME AS 

CHL PATTERN 



I 2.0 - Concluded. 

Figure 65.— Concluded. 



CONFIGURATION MEASURED 



cylinder bodies of revolution. 
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Figure 66.— Continued. 
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Figure 66.- Concluded. 
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Figure 67.- Side>foice coefncients, skfe-Toroe positions, and vapor-screen photographs for 
ogi«<e<ylinder body of revcdution with flneness-ratio-S nose. 
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for an cs«ve*cylinder body of revolution. 
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PLAN VIEW MEASURED 



Figure 69.— Effect of afterbody side strakes on side-force coefficient and side-force position for an 
ogive-cylinder body of revolution with a fineness-ratio-5 nose. 
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Figure 69.- Concluded. 
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Figure 70.— Effect of elliptic cross section and roll angle on side-force coefficient and 
side-force position (elliptic body of constant a/b). 
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Figure 70.- Concluded. 





Figure 7 1 Effect of elliptic cross section and roll angle on side-force coefficient and 
side-force position (elliptic body of variable a/b). 
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Figure 71.— Concluded. 
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(a) M^ = 0.6, Re = 4.3X10*. 

Figure 72.- Effects on side-force coefficient and side-force position of adding a wing and a 

wing plus tail to a body. 












